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Abstract 

We study the chromatic polynomials (= zero-temperature antiferromagnetic 
Potts-model partition functions) Pg(q) f° r m x n rectangular subsets of the 
square lattice, with m < 8 (free or periodic transverse boundary conditions) 
and n arbitrary (free longitudinal boundary conditions), using a transfer matrix 
in the Fortuin-Kasteleyn representation. In particular, we extract the limiting 
curves of partition- function zeros when n — > do, which arise from the crossing 
in modulus of dominant eigenvalues (Beraha-Kahane- Weiss theorem) . We also 
provide evidence that the Beraha numbers P>2, -B3, -B4, -B5 are limiting points of 
partition-function zeros as n — > 00 whenever the strip width m is > 7 (periodic 
transverse b.c.) or > 8 (free transverse b.c). Along the way, we prove that a 
noninteger Beraha number (except perhaps -Bio) cannot be a chromatic root of 
any graph. 



Key Words: Chromatic polynomial; chromatic root; antiferromagnetic Potts model; 
square lattice; transfer matrix; Fortuin-Kasteleyn representation; Temperley-Lieb al- 
gebra; Beraha-Kahane- Weiss theorem; Beraha numbers. 



1 Introduction 



The Potts model [p], |^, || plays an important role in the general theory of critical 
phenomena, especially in two dimensions ||, ^ and has applications to various 
condensed- matter systems 0. Ferromagnetic Potts models have been extensively 
studied over the last two decades, and much is known about their phase diagrams 
0, |3|] and critical exponents |5], ||, 0]. But for antiferromagnetic Potts models, many 
basic questions remain open: Is there a phase transition at finite temperature, and 
if so, of what order? What is the nature of the low-temperature phase(s)? If there 
is a critical point, what are the critical exponents and the universality classes? The 
answers to these questions are expected to be highly lattice-dependent, in sharp 
contrast to the universality typically enjoyed by ferromagnets. 

According to the Yang-Lee picture of phase transitions ||, information about 
the possible loci of phase transitions can be obtained by investigating the zeros of 
the partition function when one or more physical parameters (e.g. temperature or 
magnetic field) are allowed to take complex values. For the Potts model on a finite 
graph G, the partition function Z G (q, v) depends on the number q of Potts states and 
on the temperature-like variable v = e^ J — 1. The Fortuin-Kasteleyn representation 



|9|, |IIJ] shows that Zc(q,v) is a polynomial in q and v (see Section [2.1|) , so it makes 
sense for either or both of these variables to be made complex. In particular, the 
chromatic polynomial Pg{q) — Zg(<1, — 1) corresponds to the zero-temperature limit 
of the antiferromagnetic Potts model (J = — oo, v — — 1). 

Many investigations of the zeros of Potts partition functions in the complex q- 
and/or v -plane have been performed in the last few years, notably by Shrock and 
0, 0, 0, ID, 00, 



collaborators 
|30l |3T| 
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18, 19, 20, 21, 22, 



28, 29 



32] , |33| , [34], [35(] . The best results concern families G n of graphs for which the 



partition function can be expressed via a transfer matrix T of fixed size M x M: 



Z Gn (q,v) = tr[A{q,v)T(q,v) n ) 



M 



^2a k (q,v) X k (q, 



'1.1a} 



;i.ib) 



k=l 



where the transfer matrix T(q, v) and the boundary-condition matrix A(q, v) are 
polynomials in q and v, so that the eigenvalues of T and the amplitudes {a^} 
are algebraic functions of q and v. It then follows, using a theorem of Beraha- 
Kahane- Weiss [J7|, |SP|, [HJ, that the zeros of Zc n {q,v) accumulate along the 
curves B where T has two or more dominant eigenvalues (i.e. eigenvalues of maximum 
modulus), as well as at the isolated points where T has a single dominant eigenvalue 
\ k whose corresponding amplitude a k vanishes. See Section B72] for more details. 



For ferromagnetic Potts models on the square, triangular and hexagonal lattices, 
the exact critical curves v c (q) in the real (q, v )-plane have long been known ||. For an- 
tiferromagnetic Potts models, by contrast, there are some tantalizing conjectures con- 
cerning the critical loci, but many aspects remain obscure.Q The two best-understood 



1 For a more detailed review, see jo, Section 1]. 
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cases appear to be the square and triangular lattices: 

Square lattice. Baxter |4], 42] has determined the exact free energy (among other 
quantities) for the square-lattice Potts model on two special curves in the (q, v )-plane: 



v = ±yTq 

v = -2±y/A-q 



(1.2) 
(1.3) 



Curve (|1.2| + ) is known to correspond to the ferromagnetic critical point, and Baxter 

corresponds to the antiferromagnetic critical point, 
for q 



3 it predicts a zero-temperature 



conjectured that curve 
For q = 2 this gives the known exact value ||3| 
critical point (v c = —1), in accordance with strong analytical and numerical evidence 
~~ 1 |48|, [501 ; and for q > 3 it predicts that the putative critical point 



44, 45, M 



lies in the unphysical region (v < -1 on; complex), so that the entire physical region 
— 1 < v < lies in the disordered phase, in agreement with numerical evidence for 
q = 4 |50| . For some interesting further speculations, see Saleur |51, [52]. 

Triangular lattice. Baxter and collaborators |53], [54], |55[] have determined the 



exact free energy (among other quantities) for the triangular-lattice Potts model on 
two special curves in the (g, t>)-plane: 



v A + 3v 2 - q = 
v = -1 



;i.4) 
;i.5) 



The uppermost branch (v > 0) of curve (|1.4|) is known to correspond to the ferromag- 
netic critical point |53|, f§; and Baxter J54J initially conjectured (following a hint of 
Nienhuis ][56|) that ( |1 . 5| ) — which is the zero-temperature antiferromagnetic model, 
hence the chromatic polynomial — corresponds in the interval < q < 4 to the 
antiferromagnetic critical point. This prediction of a zero-temperature critical point 
is known to be correct for q = 2 |57|, [5£| [59| and is believed to be correct also for 

q = 4 m btl m 



On the other hand, for q = 3 this prediction contradicts the rigor- 
ous result [JS3], based on Pirogov-Sinai theory, that there is a low-temperature phase 
with long-range order and small correlation length.0 For the model (|1.5|) , Baxter |)4j 
computed three different expressions A,(g) [i = 1,2,3] that he argued correspond to 
the dominant eigenvalues of the transfer matrix in different regions V>i of the com- 
plex g-plane; in a second paper p^] he provided corrected estimates for the precise 
locations of T>i,T> 2 ,T> 3 . Unfortunately, no analogous analytic prediction is available 
for the chromatic polynomials of other two-dimensional lattices. 

One way to test the conjecture that ( |1.3| + ) is a critical curve for the square- lattice 
Potts model is to compute the partition function Z mxn (q, v) for m x n strips of 
the square lattice, investigate its zero variety in the complex (g, t>)-space, and test 
whether the zeros of Z mxn (q, v) appear to be converging to (101+) as m,n^ oo. Here 
we shall carry out this program for the zero-temperature antiferromagnetic model 



2 A Monte Carlo study of the q = 3 model found strong evidence for a first-order transition to 
an ordered phase at @J sa —1.594 p3. 
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{y = — Using a transfer matrix in the Fortuin-Kasteleyn representation [66 



we 



shall compute the chromatic polynomials P mX n((l) for m x n square-lattice strips of 
width m < 8 (free or periodic transverse boundary conditions) and arbitrary length 
n (free longitudinal boundary conditions). In particular, we shall extract the limiting 
curves B of partition-function zeros when n —>■ oo, which arise from the crossing 
in modulus of dominant eigenvalues in accordance with the Beraha-Kahane-Weiss 
theorem.^ Finally, we shall attempt to understand the behavior of these limiting 
curves as m — > oo. Of course, there is little doubt in this case that they will converge 
to the critical point of the zero-temperature model, q c = 3; but it is illuminating to see 
this convergence explicitly and to view the critical point q c = 3 as simply one (real) 
point on a complex critical curve. Not surprisingly, we find for this critical curve 



a shape that is qualitatively similar to that found by Baxter [||IJ for the triangular 
lattice, with the zero-temperature critical point lying now at q c = 3 rather than 
q c = 4. 

A special role in the theory of chromatic polynomials appears to be played by 
the Beraha numbers B n = 4cos 2 (-7r/n) [see Table [TJ for the first few B n }. As we shall 
show in Section |2.3| , a noninteger Beraha number (except possibly Bio) cannot be a 
chromatic root of any graph. Nevertheless, Beraha |J7| observed that planar graphs 
frequently have chromatic roots very near one or more of the B n .^ Indeed, Beraha, 



Kahane and Weiss |38], found families of planar graphs that have chromatic roots 
converging to B 5 , B 7 or B w ^\ 

Here we shall provide additional curious evidence in favor of the idea that chro- 
matic roots tend to accumulate at the Beraha numbers. We find empirically (at least 
for m < 8) that on a square-lattice strip of width m with either free or periodic 
transverse b.c, there is at least one vanishing amplitude ati(q) at each of the first 
m Beraha numbers B2, . . . , B m+ i (but not higher ones). Assuming that this behav- 
ior persists for all m, in the limit m — > 00 all the Beraha numbers will be zeros of 
some amplitude. Moreover, in all the cases except m = 7, 8 with free transverse b.c. 
(where our computer power gave out) and m = 8 with periodic transverse b.c. (see 

3 In future work |65| we plan to extend this analysis to the antiferromagnetic model at nonzero 
temperature (—1 < v < 0). See also [ p0| . 

4 Here we follow in the footsteps of Shrock and collaborators |l7], [ll| [32j , who have been carrying 
out this program using a generating-function approach that is equivalent to transfer matrices; they 
determine the recurrence relations by repeated use of the deletion-contraction identity. In particular, 
Shrock et al. have computed the transfer matrices for square-lattice strips of width m < 5f and 
m < 6p (leading to matrices of size up to 7 x 7), and have computed the limiting curves B of 
partition-function zeros for m < 4f and m < 5p. [Here the subscript F (resp. P) denotes free 
(resp. periodic) boundary conditions.] By explicit use of transfer matrices, we are able to automate 
the former calculation and ha ndle m uch larger transfer matrices (here up to 127 x 127); and using 
the resultant method (Section |4.1.l| ) we are able to detect small gaps and other fine details in the 
limiting curves. 



5 For B5 this was observed earlier by Berman and Tutte |68 . 

6 The graphs in question are 4p x tif, 5p x n-p and 2f x np strips of the triangular lattice (with an 
extra vertex adjoined at top and bottom in the first two cases). In the case of B§ and Br, Beraha, 
Kahane and Weiss |39f| proved that there are even real chromatic roots converging to them. 
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Section for a discussion on this point), we verified that the vanishing amplitude 
corresponds to the eigenvalue obtained by analytic continuation in q from the one 
that is dominant at small real q (e.g. at q = 1), in agreement with a conjecture of 
Baxter J55, p. 5255]. Thus, the first few Beraha numbers — namely, those (up to at 
most B m+ i) that lie below the point qo(m) where the dominant-eigenvalue-crossing 
locus B intersects the real axis — correspond to the vanishing of a dominant ampli- 
tude and hence (via the Beraha-Kahane-Weiss theorem) to a limit point of chromatic 
roots, while the remaining Beraha numbers do not. As the strip width m grows, this 
crossing point qo(m) increases and presumably tends to a limiting value go(°°); f° r 
the square lattice, we expect qo{oo) to lie somewhere around 2.9, i.e. strictly between 
B 5 and B G . Therefore, for all sufficiently large strip widths, we expect the Beraha 
numbers B 2 , B 3 , B4, B 5 — but not higher ones — to be limiting points of chromatic 
roots. Our data confirm (at least up to m = 8) that B 2 ,B 3 , B A are limiting points of 
zeros for all widths m > 4, and that B 5 is a limiting point of zeros for all widths m > 7 
(cylindrical b.c.) or in > 8 (free b.c). This scenario for the accumulation of chromatic 
roots at some of the Beraha numbers was set forth by Baxter [55 and elaborated 



by Saleur |5T|. For further speculations on the special role of the Beraha numbers in 
the Potts model, and especially for the chromatic polynomials of planar graphs, see 

H, §9], 0, o HI 0, 0, SiiSS BUI- 

The plan of this paper is as follows: In Section § we review the Fortuin-Kasteleyn 
representation, the Beraha-Kahane-Weiss theorem, and some algebraic number the- 
ory related to the Beraha numbers. In Section ^| we explain how to construct transfer 
matrices for the Potts model in the spin representation and in the Fortuin-Kasteleyn 
representation, and we compute the dimensions of these transfer matrices. In Sec- 
tion [| we discuss the general properties of the dominant-eigenvalue-crossing curves B 
and the isolated limiting points of zeros. In Sections [5] and |6| we present our numeri- 
cal results for square-lattice strips with free and cylindrical boundary conditions. In 
Section [7] we analyze the theoretical import of our calculations, and discuss prospects 
for future work |86|, |8 



8,B5 



2 Preliminaries 



In Sections EO and |2.2| we review some well-known facts about the Fortuin- 
Kasteleyn representation and the Beraha-Kahane-Weiss theorem, which will play 
a fundamental role in the remainder of the paper. We also use these sections to set 



the notation. In Section 2.3 we discuss some algebraic number theory related to the 



Beraha numbers; this section contains a few new results, notably Corollary 2.4. 



2.1 Fortuin-Kasteleyn representation of the Potts model 

Let G = (V, E) be a finite undirected graph with vertex set V and edge set E, let 
{J e }eeE be a set of couplings, and let q be a positive integer. Then the q-state Potts 
model on G with couplings {J e } is, by definition, the canonical ensemble at inverse 
temperature j3 for a model of spins {o~ x } x£ v taking values in the set {1,2, ... ,q}, 
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interacting via a Hamiltonian 



#(W) = - X) J AWv) (2-1) 

e={xy)eE 

where S is the Kronecker delta. The partition function is thus 

Z G (q,{v e }) = X II + (2.2) 

{<r} e=(xi/>g£; 



where we have written 



e p ' h - 1 . (2.3) 



A coupling J e (or v e ) is called ferromagnetic if J e > (t> e > 0) and antiferromagnetic 
if — oo < J e < (— 1 < v e < 0). The g-coloring problem, in which adjacent spins 
are required to take different values, corresponds to the zero-temperature limit of the 
antiferromagnetic Potts model (namely J e = — oo, v e = —1). 

In fact, Zc{q, {v e }) is the restriction to positive integers q of a polynomial in q and 
{v e } (with coefficients that are in fact or 1). To see this, expand out the product 
over e G E in (|2.2p, and let E' C Ebe the set of edges for which the term v e S(a x , a y ) is 
taken. Now perform the sum over configurations {a}: in each connected component 
of the subgraph (V, E') the spin value a x must be constant, and there are no other 
constraints. Therefore, 

Z G (q,{ve}) = E^II^' ( 2 - 4 ) 

E'CE e£E' 

where k(E') is the number of connected components (including isolated vertices) 
in the subgraph (V,E'). The expansion ( |2.4|) was discovered by Birkhoff |89| and 
Whitney |Kj for the special case v e — — 1 (see also Tutte 0, [9^]); in its general 
form it is due to Fortuin and Kasteleyn M, ITH] (see also pi). We shall henceforth 



take ( |2.4| ) as the definition of Zc(q, {v e }) for arbitrary complex numbers q and {f e }. 
When v e takes the same value v for all edges e, we write Zc{q,v). When v e = — 1 
for all edges e, this defines the chromatic polynomial Pc{q)- Note that the chromatic 
polynomial Pc{q) of any loopless graph G is a monic polynomial in g with integer 
coefficients.^ See |94| , |95| for excellent reviews on chromatic polynomials, and |K| for 
an extensive bibliography. 

2.2 Beraha— Kahane— Weiss theorem 

A central role in our work is played by a theorem on analytic functions due to 
Beraha, Kahane and Weiss [B5|, [J7], and generalized slightly by one of us pi 



The situation is as follows: Let D be a domain (connected open set) in the complex 

7 A loop, in graph-theoretic terminology, is an edge connecting a vertex to itself. Obviously, if G 
has a loop, then its chromatic polynomial is identically zero. 
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plane, and let ai, . . . , Ai, . . . , Am (M > 2) be analytic functions on D, none of 
which is identically zero. For each integer n > 0, define 



M 



f n (z) = Y,<*k{z)\ k {z) n . (2.5) 

k=l 

We are interested in the zero sets 

Z(f n ) = {zeD: f n (z) = 0} (2.6) 

and in particular in their limit sets as n — > oo: 

liminf Z(f n ) = {z G D: every neighborhood U 3 z has a nonempty intersection 

with all but finitely many of the sets Z(f n )} (2-7) 

limsupi?(/ n ) = {z G D: every neighborhood U 3 z has a nonempty intersection 

with infinitely many of the sets Z(f n )} (2-8) 

Let us call an index k dominant at z if |Afc(z)| > |A;(z)| for alii (1 < I < M); and let 
us write 

-D/c = {z E D: k is dominant at 2} . (2-9) 
Then the limiting zero sets can be completely characterized as follows: 



Theorem 2.1 [|3q , |37| , |38| , |39| , Lei D be a domain in C, and letai,. . . , olm, Ai, . . . , Am 

> 2) be analytic functions on D, none of which is identically zero. Let us further 
assume a "no- degenerate- dominance" condition: there do not exist indices k 7^ k! 
such that Afc = ouXjc' for some constant u with \u\ = 1 and such that Dk (= Dk> ) has 
nonempty interior. For each integer n > 0, define f n by 

M 

f n {z) = Y^ a ^z)X k {z) n . 
k=i 

Then lim inf Z(f n ) = limsup Z(f n ), and a point z lies in this set if and only if either 

(a) There is a unique dominant index k at z, and a k {z) = 0; or 

(b) There are two or more dominant indices at z. 

Note that case (a) consists of isolated points in D, while case (b) consists of curves 
(plus possibly isolated points where all the A^ vanish simultaneously). Beraha- 
Kahane- Weiss considered the special case of Theorem 2~1 in which the /„ are polyno- 
mials satisfying a linear finite-order recurrence relation, and they assumed a slightly 
stronger nondegeneracy condition. (This is all we really need in this paper.) Hence- 
forth we shall denote by B the locus of points satisfying condition (b). 

In the applications considered in this paper, the functions f n will be of the form 
f n (z) = tr[A(z)T (z) n ] where the M x M matrices T(z) and A(z) are polynomials in 
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z; therefore, the functions Xk(z) [which are the eigenvalues of T(z)] and aik(z) will be 
algebraic functions of z, i.e. locally analytic except at isolated branch points. So one 
can cover the complex plane minus branch points by a family of simply connected 
domains D, and then apply Theorem [2.1| separately to each such domain D. The 
branch points are not covered by this analysis, but they will always be endpoints of 
curves of type (b), hence also limit points of zeros. 

It is interesting to ask about the rate at which the zeros of /„ converge to the 
limit set. In case (a), it is easy to see that the convergence is exponentially fast. In 
case (b), simple expansions suggest that the rate is 1/n near regular points of the 
curve B, and 1/n 2 at branch points (see Sections [4.2. 1| and |4.2.2|) . 

In checking for isolated limiting points of zeros [case (a)], the following result will 
be useful (see Section 14.31): 



M 



Lemma 2.2 ||39|| Suppose that f n = ^2 a k A£, and define 



k=l 



D 



( fo fl'" /m-i \ 
fl f'2 ■ ■ ■ /m 

\ /m-i /m • • • /2M-2 / 



(2.10) 



Then 



M 



det D = ]Ja k (Xj - Xi) 2 . 

k=l l<i<j<M 



(2.11) 



Proof. It is not difficult to see that D = A T diag(«i, . . . , «m)A where A is the 
M x M Vandermonde matrix A^- = A| - . Lemma |2.2| then follows from the well- 
known formula 

det A = I [ [Xj - Xi) (2.12) 

l<i<j<M 

for the Vandermonde determinant. ■ 

We remark that the product n i<3 -(Aj — Xi) 2 is the discriminant of the characteristic 
polynomial of the transfer matrix: see ( |4.2j ) below. 



2.3 Beraha numbers 

We recall [0] that a complex number ( is called an algebraic number (resp. an 
algebraic integer) if it is a root of some monic polynomial with rational (resp. inte- 
ger) coefficients. Corresponding to any algebraic number (, there is a unique monic 
polynomial p with rational coefficients, called the minimal polynomial of ( (over the 
rationals), with the property that p divides every polynomial with rational coefficients 
having ( as a root. (The minimal polynomial of ( has integer coefficients if and only 
if ( is an algebraic integer.) Two algebraic numbers ( and (' are called conjugate 



s 



if they have the same minimal polynomial. Conjugacy is obviously an equivalence 
relation, and so divides the set of algebraic numbers into equivalence classes. 

Examples. 1. The numbers B 5 = (3+Vb)/2 and B* = (3-y/E)/2 are conjugates, 
as they have the same minimal polynomial p(x) = x 2 — 3x + 1. 

2. The numbers 2 1//3 , 2 1 / 3 e 27 ™ //3 and 2 1/ ' 3 e 47n/ ' 3 are all conjugates, as they have the 
same minimal polynomial p(x) = x 3 — 2. 

3. If C is a primitive n th root of unity, the minimal polynomial of ( is the cyclotomic 
polynomial 

$n(x) = J] (x - e 2 ^ n ) , (2.13) 

1 < k < n 
gcd(fc, n) = 1 

where the product runs over all positive integers k < n that are relatively prime to n 
(see e.g. 0, Section 3.7] for a proof). In particular, $ n is a monic polynomial with 
integer coefficients, irreducible over the field of rational numbers, of degree 

deg$ n = (p(n) (2.14) 

where ip(n) is the Euler totient function (i.e. the number of positive integers k < n 
that are relatively prime to n). Thus, all the primitive n th roots of unity are mutually 
conjugate. 

4. For n > 2, let us define the generalized Beraha numbers 

flW = 4 cos 2 — = 2 + 2 cos , 2.15 

n n 

and let us call a primitive n th generalized Beraha number in case k is relatively 
prime to n. Then the minimal polynomial of B n (or of any primitive Bn ) isQ 

Pn(x) = I! ( X - B n ] ) ■ ( 2 - 16 ) 

1 < k < n/2 
gcd(fc, n) = 1 

In Table [l] we show the p n for 2 < n < 16. p n is a monic polynomial with integer 
coefficients, irreducible over the field of rational numbers, of degree 

_ f 1 for n = 2 , n , 

degP " = W(n)/2 forn>3 (2 ' 17) 



The proof of ( 2.16 ) requires some elementary Galois theory |98|, Chapter 3] J99[] : Fix n > 3 and 
C = e 27ri /™ ; an( j \ e i p \, e iy^q extension field Q(C)- The irreducibility of the cyclotomic polynomial 
<&„ implies that F has dimension <p(n) over Q, and that its Galois group G = Gal(F/Q) is the 
abelian group {ck}i<k<n,ecd(k,n)=i where <7fc(r) = r for r € Q and <Tfc(C) = C fc - Now define c = 
2cos(27r/n) = ( + ( £ R. By repeated application of the equation ( 2 — c( + 1 = it follows that 
Q(C) = Q( c ) © CQ( c )j so that F has dimension 2 over Q(c). Therefore, Q(c) has dimension ip(n)/2 
over Q, so that the minimal polynomial P(x) of c over Q has degree ip(n)/2. Now, for every A; that 
is relatively prime to n, let us define Ck = crfe(c) = ^ fe + ( k ; there are clearly ip(n)/2 distinct such 
numbers Cfc (note that c„_fe = c/.). And we have P(ck) = P{cfk{c)) = o"/-(P(c)) = since is a field 
automorphism and P has rational coefficients. It follows that P{x) = Ili<fe<n/2 gcd(fc n)=i( :E — c *;)- 
A trivial shift x — > a; — 2 gives the same result for the £?„ ; = 2 + c&. We thank Dan Segal for 
explaining this proof to us. 
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Thus, all the primitive n generalized Beraha numbers are mutually conjugate. 

The point of this digression into algebraic number theory is, of course, that the 
chromatic polynomial Pg(q) of any loopless graph G is a monic polynomial with 
integer coefficients. Therefore, all the chromatic roots of G are algebraic integers; 
and if ( is a chromatic root of G, then so are all its conjugates. We therefore have: 

Proposition 2.3 (a) Suppose that ( has a conjugate lying in one of the intervals 
(—oo,0) ; (0,1) or (1,32/27]. Then ( is not a chromatic root of any loopless 
graph. 

(b) Suppose that ( has a conjugate lying in the interval [5, oo). Then £ is not a 
chromatic root of any loopless planar graph. 

(c) Suppose that ( has a conjugate lying in the interval (1,2). Then ( is not a 
chromatic root of any plane near-triangulation. 

(d) Suppose that ( has a conjugate lying in the interval (2,2.546602...), where 
2.546602 . . . is shorthand for the unique real solution of q 3 — 9q 2 + 29q — 32 = 
(which is the nontrivial chromatic root of the octahedron). Then ( is not a 
chromatic root of any plane triangulation. 



Proof. This follows immediately from the facts that there are no chromatic roots 
in the relevant intervals: see ref. |95] for (— oo,0) and (0,1), ref. |100|| for (1,32/27], 



refs. [fTTH], [Iig, |TU§j for [5, oo), refs. JTDTj, TO PU5| for (1,2), and refs. |TU6J for 



(2,2.546602 



Examples. 1. A special case of Proposition |2.3j (a) is the well-known fact |55 



p. 23] that £> 5 = (3 + y/E)/2 is not a chromatic root of any (loopless) graph, since its 
conjugate = (3 — v / 5)/2 lies in (0, 1) [their minimal polynomial is q 2 — 3q + 1] . We 
shall prove a vastly stronger result in Corollary |2.4| below. 

2. The number 1 + \/2 is not a chromatic root of any (loopless) graph, since its 
conjugate 1 — a/2 lies in (— oo, 0) [their minimal polynomial is q 2 — 2q — 1]. 

3. For n > 5, none of the numbers 2 1//n e 27rfci//n (k integer) is a chromatic root 
of any (loopless) graph, since their conjugate 2 1//ra lies in (1,32/27] (their minimal 
polynomial is q n — 2).0 

4. The numbers 2 + \j 2\[2 — 2 and 2 ±iy^2\^2 + 2 are not chromatic roots of any 
(loopless) graph, since their conjugate 2- A/2A/2-2 lies in (1,32/27] (their minimal 
polynomial is q A - 8q 3 + 28q 2 - 48g + 28). 

We now state and prove the promised generalization of the fact that B 5 is not a 
chromatic root: 



9 See also [102, 103 for related results on zero- free intervals. 

10 This barely fails for n = 4. It is amusing to note that 32/27 as 2 1 / 4 is equivalent to (3/2) 12 
which underlies the construction of the well-tempered scale. 
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Corollary 2.4 A noninteger primitive n th generalized Beraha number Bn is not a 
chromatic root of any graph, except possibly when n — 10. 

This is an immediate consequence of the case (0, 1) of Proposition pH5|(a) together 
with the following number-theoretic lemma: 

Lemma 2.5 Let n = 5, 7, 8, 9 or n > 11. Then there exists an integer k, relatively 
prime to n, in the interval n/3 < k < n/2. 



Proof. If n is a prime > 5, the result is trivial (take k = \n/2\). The cases n — 8, 12 
can be verified by hand (take k = 3,5, respectively). Otherwise, note first that it 
suffices to find k in the interval n/3 < k < 2n/3 (if k > n/2, then replace k by 
n — k). So let p be the least prime divisor of n; define r = \p/3] > 1, so that 
rn/p > n/3 > (r — l)n/p; and let k be either rn/p + 1 or rn/p + 2. Now the only 
possible common prime factor of k and n is p: for if p' ^ p were to divide n, then it 
would also divide n/p; but if p' also divided k, then it would divide k — rn/p = 1 or 
2, which is impossible since p' > 3. But p cannot divide both choices of k. Therefore, 
this construction allows k to be chosen relatively prime to n; and k lies in the interval 
(n/3, 2n/3) provided that n > 12 (when p = 2), n > 6 (when p = 3) or n > 15 (when 
p>5). I 

Remarks. 1. As mentioned earlier, the case n = 5 of Corollary |2.4| is ancient 



folklore; the case n = 7 was stated by Tutte ]72], p. 372]. Beraha, Kahane and Weiss 
|39| , p. 53] asserted that the argument for B 5 "can be extended without much difficulty 
to the case of arbitrary nonintegral B n " , apparently overlooking the problem with Biq. 

2. The exceptional case n = 10 is very curious. We do not know whether B w = 
(5 + y/E)/2 w 3.6180339887 and B* = (5 - v / 5)/2 « 1.3819660113 can be chromatic 



roots. But Proposition |2.3| (c) shows that -Bio is not a chromatic root of any plane 
near-triangulation . Note also that when G is a plane triangulation having n vertices, 
Tutte's "golden identity" |70J |§, pp. 26-27] 



Pg{B w ) = (r + 2)r 3n - 10 P G (B 5 ) 2 (2.18) 

where r = B 5 — 1 is the golden ratio, together with the fact that P G (B 5 ) ^ 0, yields 
the slightly stronger result Pc(B w ) > 0. 



3 Transfer- Matrix Theory 

In this section we briefly review the transfer-matrix formalism for the Potts model, 
both in the spin representation (Section |3.1| ) and in the Fortuin-Kasteleyn represen- 
tation (Section |3.2|) . Thereafter we use only the Fortuin-Kasteleyn representation, 
which allows us to perform computations for noninteger q. We conclude by comput- 
ing the size of the FK-representation transfer matrix for the square and triangular 
lattices with free longitudinal boundary conditions and either free or periodic trans- 
verse boundary conditions (Section |3~3|). 
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3.1 Transfer matrix in the spin representation 

Consider a graph G n = (V n , E n ) consisting of n identical "layers" , with connections 
between adjacent layers repeated in a regular fashion. Then the partition function of 
the Potts model on G can be written in terms of a transfer matrix. 

To make this precise, let V° be the set of vertices in a single layer, let E° be the 
set of edges within a single layer (we call these horizontal edges), and let E* be the 
set of edges connecting each layer to the next one (we call these vertical edges) . Note 
that E° is a set of unordered pairs of elements of V°, while E* is a set of ordered 
pairs of elements of V°. The vertex set of the graph G n is then 

V n = V° x {1, . . . , n} = {(x, i): x G V° and 1 < i < n} , (3.1) 

while the edge set is either 

£free = ^oriz I I "y 1 £vert ^) 

i=l i=l 

for free longitudinal boundary conditions or 

El er = U E* OTiz I J U EJ crt (3.3) 
i=i ^-^ i=i 

for periodic longitudinal boundary conditions, where 

^ hOTiz = {{(x,i),(x',i)): (xx')EE } (3.4) 

and 

EJ crt = {((x,i),tf,i + l)): (x, x') G E*} (3.5) 

and of course layer n+ 1 is identified with layer 1. We also assume that the couplings 
are identical from layer to layer: that is, we are given weights {v e } ee E°uE*, and we 
define the edge weights for G n by 

u<(*,i),(x',i)> = v {xx , } for (xx') E E° (3.6a) 
V(( x ,i),( x >,i+i)) = v {XjX) for (x, x 1 ) G E* (3.6b) 

We now fix an integer q > 1, and let E = {l,...,q} v ° be the space of spin 
configurations on a single layer; we denote a generic such configuration by a = 
{°~x}xev°- We then define the matrices H and V, which encode the Boltzmann factors 
corresponding to horizontal and vertical edges, respectively: 

H(<t',<t) = 5{(t,(t') Y[ [l + v {xx/) 5(a x ,o- x ,)} (3.7) 

(xx')eEO 

V(<7» = n [l+v {x>xf) 5(cr x ,a' x/ )] (3.8) 

(x,x')€E* 
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The partition function of the Potts model on G n is then 

Z G ^(q,{v e }) = l T H(VH) n ~ 1 l (3.9) 

= tr[(VH)«] (3.10) 

where 1 is the vector whose entries all equal 1, and T denotes transpose. The transfer 
matrix is thus 

T = VH . (3.11) 

It has size q m x q m , where m — \V°\ is the number of sites in each layer. 

The horizontal matrix H is diagonal in the spin basis, hence sparse; its computation 
takes a time of order \E°\q m < m 2 q m . But since the vertical matrix V is dense, its 
computation takes in general a time of order \E*\q 2m . This is a severe constraint on 
the practical applicability of the transfer-matrix method. However, when the graphs 
G n are planar, V can be written as a product of sparse matrices that correspond 
to the replacement of one site on layer i by the corresponding site on layer i + 1, 
so that its computation takes only a time of order mq m . (Indeed, this can be done 
for an arbitrary planar graph, whether or not it consists of repeated layers ||79|| . ) 
This factorization of V can also be performed for some non-planar graphs, including 
all those in which the single-layer vertex set V° can be ordered in such a way that 
(x, x') G E* implies x >z x'-Q (This latter situation includes, in particular, all the 
lattices Z d , since the only vertical edges are of the form (x, x). On the other hand, it 
excludes the triangular lattice with periodic transverse boundary conditions, which 
needs a separate treatment.) For simplicity we describe this construction only for the 
cases of greatest practical interest, namely when G n is a square or triangular lattice 
with free or periodic transverse boundary conditions. The single-layer vertex set is 
thus V° = {1, . . . , m); the horizontal edge set is either 

E^ ee = {(x,x + l): l<x<m-l} (3.12) 

for free transverse boundary conditions or 

E° er = {(x,x + l): l<x<m} (3.13) 

for periodic transverse boundary conditions (where site m + 1 is of course identified 
with site 1); and the vertical edge set is either 

E° SQ = {(x,x): l<x<m} (3.14) 

for the square lattice, 

^TRI.free = X Y 1 < X < m } IJ { ( X + 1 > X ) ■ I < X < TU - 1} (3.15) 

for the triangular lattice with free transverse boundary conditions, or 

^TRi.per = {( x , x )- 1 < x < m} |J { (x + 1 , x) : 1 < x < m} (3.16) 

11 This is equivalent to requiring that the directed graph (V°,E**) be acyclic, where E** = 
{(x,x') € E*: x ^ x'}. 
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for the triangular lattice with periodic transverse boundary conditions. Note that the 
diagonal edges in (|3.15|) - (|3.16|) point "southeast-northwest". 



We now define q m x q m matrices and S xx i by 



D x {a',a) = \[8{a v ,a' y ) (3.17a) 
} xx r(cr',cr) = 8(a x ,a x/ )8(cr,cr') (3.17b) 



Thus, D x (for "detach" or "disconnect") is of the form I ® I ® ■ ■ ■ ® E ® ■ ■ ■ ® I <g> I, 
where I is the q x q identity matrix, E is the q x q matrix whose entries are all 1, 
and the E factor occurs at position x; informally, D x disconnects the two rows at 
site x. Since I and E commute, all the operators D x commute among themselves. 
Likewise, S xx i identifies ("joins") spins x and x' in a single layer. Since the operators 
S xx i are diagonal in the spin basis, they also commute among themselves. Finally, 
commutes with S x i x ii whenever x is different from both x' and x" . 
We now define matrices 

Pi = + D, (3.18) 

Qj = I + v {jjj+1) J jd+1 (3.19) 

R j = 1 + ( 3 - 20 ) 

corresponding to the Boltzmann factors for vertical, horizontal and diagonal edges, 
respectively.^ The horizontal and vertical parts of the transfer matrix are then 

Hsq c = Htri = Q™- 1 ' ' ' Q2Q1 (3.21a) 

Hsq = QmQm-i • • • Q2Q1 (3.21b) 

and 

Vsq 6 = Vgq = P m P m -i • ■ • P2P1 (3.22a) 

^TRI = Pm^m— lPm— 1 ' ' ' R2P2R1P1 (3.22b) 

where the action of these matrices should always be read from right to left. Note that 
Rj corresponds to a diagonal edge whenever it is applied after Pj and before Pj+i- 

For the triangular lattice with periodic transverse boundary conditions, a slight 
trick is needed in order to treat correctly the last diagonal bond joining columns 
m and 1 |55|. The idea is to work with q m+l x q m+1 matrices indexed by spins 
<t = (<7i, . . . , <7 m+ i), and include in the horizontal matrix an operator that identifies 
°~ m +i with <Ti on each row. Thus, 

^TRI = Jm+l,lQmQm-l * " ' Q2Q1 (3.23) 
^TRI = Pm+lPmPmRm-lPm-1 " " " R2P2R1P1 (3.24) 



12 In this sentence we are of course using the term "vertical edges" in its ordinary geometrical 
sense, not in the technical sense defined previously. We follow the notation of Baxter |B5|, Section 2] 
with some modifications. 
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In particular, the vertical matrix V^rt is just V^rt with m replaced by m + 1. 

Remarks. 1. The rewriting (3.21) of the horizontal matrix H is of course unnec- 
essary, as H was already diagonal (hence sparse) in the spin basis. But (3.21) will be 
useful when we use the partition basis (see below). 

2. The order of operators in (3.21), ( |3.22aj ) and ( |3.23j ) is irrelevant, as all the 
operators in question commute. Only in (3.22b) and (|3.24|) is the operator ordering 
crucial. 

3. The operators ei, . . . , e 2m _i defined by 

e 2j _ 1 = q~ 1/2 Dj for 1 < j < m (3.25a) 
e 2j = q 1/2 l jJ+1 for 1 < j < m - 1 (3.25b) 



satisfy the Temperley-Lieb algebra ||107| , £|, [79| 

e\ = q 1/2 e i (3.26a) 

eie i±1 ei = (3.26b) 

titj = tjti for \i — j\ > 1 (3.26c) 



All these operators act, of course, in the g m -dimensional vector space consisting of 
all functions of <7i, . . . , cr m .Q But for the case of free longitudinal boundary conditions 
(p.9|), we need only consider the subspace spanned by those vectors that can be 
obtained from the constant vector 1 by action of the operators H and V. All such 
functions are sums of products of delta functions 5(a x , a y ), and we can take as a basis 
the functions vp associated to a partition V = {Pi, . . . , Pk} of {1, . . . , m} by 

vp(cri, . . . , a m ) ] [ S(a x ,a y ). (3.27) 

PeV x, y eP 

For example, for m = 3 we have the five basis vectors 

V{{i},{ 2 },{3}} = 1 (3.28a) 

V{{i, 2 },{3}} = S(ai,cr 2 ) (3.28b) 

V{{i, 3 },{2}} = S(a u cr 3 ) (3.28c) 

V{{!},{ 2 ,3}} = 5(a 2 ,a 3 ) (3.28d) 

V{ {1,2,3}} = 5{(Ti,a 2 ,a 3 ) (3.28e) 

Indeed, when the graph G is planar, it is not hard to see on topological grounds 
that only non-crossing partitions can arise. (A partition is said to be non-crossing 
if a < b < c < d with a, c in the same block and b, d in the same block imply that 



13 Or, in the case of the triangular lattice with periodic transverse b.c, in the g m+1 -dimensional 
space of all functions of a\, . . . , cr m +i- 
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a, b, c, d are all in the same block.) Moreover, spatial symmetries can further restrict 
the subspace. Finally, when the horizontal couplings v/ xx >\ are all equal to — 1 - 
which is the case for the chromatic polynomial — then the horizontal operator H is 
a projection, and we can work in its image subspace by using the modified transfer 
matrix T' = HVH in place of T = VH, and using the basis vectors 

w P = Hvp (3.29) 

in place of v-p. Note that wp = if V has any pair of nearest neighbors in the same 



block. In Section |3.3| we shall compute the dimensions of all these subspaces. 

The action of the operators and J xx > on the basis vectors v-p is quite simple. 
As one might expect, S xx i joins sites x and x', i.e. 

Jzi'V-p = y-puxx 1 (3.30) 

where V»xx' is the partition obtained from V by amalgamating the blocks containing 
x and x' (if they were not already in the same block). detaches site x from the 
block it currently belongs to, multiplying by a factor q if x is currently a singleton: 

\ qv v if M 6 V 

where V \ x is the partition obtained from V by detaching x from its block (and thus 
making it a singleton). 

At the final stage we need to compute the inner products l T vp, which are easy: 

l T Vp = (3.32) 

where \V\ is the number of blocks in V . 



To summarize: The Potts-model partition function ( |3.9j ) or ( |3.10| ) can always be 



computed by the transfer-matrix method working in the spin basis /(<7i, . . . , <7 m ), 
which has dimension q m . Of course, each value of q has to be treated separately; and 
it goes without saying that q must be a positive integer. On the other hand, for free 
longitudinal boundary conditions, there is an alternative method for computing the 
partition function ( |3.9| ), using the partition basis v-p. In this formalism q appears 
only as a parameter [in (|3.31| ) and ( [3.321) 1; and since the final result will obviously 



be a polynomial in q, it must perforce be equal to the Fortuin-Kasteleyn partition 
function ( |2.4|) . This latter formalism is thus essentially equivalent to constructing 
the transfer matrix directly in the Fortuin-Kasteleyn representation, as we shall see 
explicitly in the next subsection. 

3.2 Transfer matrix in the Fortuin-Kasteleyn representation 

Let us now temporarily forget the Potts spin model — and in particular forget the 
transfer-matrix formalism constructed in the preceding subsection — and try instead 
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to devise a transfer-matrix method for computing the Fortuin-Kasteleyn partition 
function ( |2.4|) when the graph G has a layered structure ( |3.1|) -(3.6). What makes this 
a bit tricky is, of course, the nonlocality of the factor q k ( E > in ( |2.4| ). The technique 
for handling this nonlocality was foreshadowed in the work of Lieb and Beyer ||108| 



on percolation and was made explicit (for the case of the chromatic polynomial) in 
the work of Biggs and collaborators [ |109| , |110 , 11 1|. In the physics literature, this 



approach was first used (to our knowledge) by Derrida and Vannimenus ||112|| in their 
study of percolation, and was applied to the g-state Potts model (and explained very 
clearly) by Blote and Nightingale ; it was subsequently employed by several groups 
||113| , |114| , |115| , |116|| . Here we limit ourselves to giving a brief summary. 



The basic idea is to build up the subgraph E' C E layer by layer. At the end we 
will need to know the number of connected components in this subgraph; in order 
to be able to compute this, we shall keep track, as we go along, of which sites in 
the current "top" layer are connected to which other sites in that layer by a path of 
occupied bonds (i.e. bonds of E') in lower layers. Thus, we shall work in the basis of 
connectivities of the top layer, whose basis elements \-p are indexed by partitions V 
of the single-layer vertex set V°. [The reader is reminded to forget ( |3.27| ) for the time 



being. The v-p should here be thought of simply as abstract basis elements indexed 
by partitions V.] The elementary operators we shall need are the join operators 

Jzz' v -p = v Ptix' (3.33) 
(note that all these operators commute) and the detach operators 

\qv P ii{x}EV 

where V»xx' and V\x were defined previously (don't forget those definitions!). The 
horizontal transfer matrix is then 

H = ] J [1 + V( XX >) Sxx>) ■ (3.35) 

(xx')£E° 

The vertical transfer matrix is slightly more complicated: 

Vv„ = J2 ^ ( II J v v\e (3-36) 

ECE* \(x,x')€E ) 

where A(V,E) is the number of "abandoned clusters", i.e. the number of blocks 
P G V such that no vertex in P is an endpoint of an edge in E\ and V\E is the 
partition of V° in which vertices x', y' lie in the same block if and only if there exist 
vertices x, y in the same block of V such that both (x, x') and (y, y') lie in E. 

In many cases (including all planar graphs), V can be written as a product of 
sparse matrices that correspond to the replacement of one site on layer i by the 
corresponding site on layer i + 1; and these sparse matrices have a simple expression 
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in terms of join and detach operators. Suppose, for concreteness, that the single-layer 
vertex set V° can be ordered in such a way that (x, x') G E* implies x >z using 
this ordering, let us number the sites of V° as 1, . . . , m. Then 

V = P m Rm-l,mPm-l I [ [ ^m-2,y J Pm-2 ' ' ' I \\ ^ 2 ' y J ^ 2 ( IT ) ^ 1 ' 

\y>m-2 J \y>2 J \y>l J 

(3.37) 

where 

P x = v (x , x) I + D x (3.38) 

^x,y = I ~\~ V{y,x)^xy (3.39) 

(The triangular lattice with periodic transverse b.c. is handled using the trick dis- 
cussed in the preceding subsection.) 

Finally, the partition function for free longitudinal boundary conditions can be 
obtained by sandwiching the transfer matrix between suitable vectors on right and 
left: 

Z G *ee{q,{v e }) = u^VH^-Vid , (3.40) 

where "id" denotes the partition in which each site x G V° is a singleton, and u T is 
defined by 

u T v P = q m . (3.41) 

Of course, it will not have escaped the reader's notice that all the formulae in 
this subsection are identical to those developed in the preceding subsection; only 
the interpretation is different. In particular, the operators D x and J xx i in the FK 
representation act in precisely the same way as the operators of the same name act 
in the spin representation with respect to the partition basis. 

Let us observe, in conclusion, that the Potts model with/ree longitudinal boundary 
conditions is much easier to handle than that with periodic longitudinal boundary 
conditions. We can understand the difficulty posed by periodic longitudinal b.c. from 
two complementary points of view. In the spin representation, the trouble is that, 
because of the trace in (|3.10| ), we must compute the action of T in the entire q m - 
dimensional space of functions of 01, . . . , c m > not merely in the subspace spanned 
by sums of products of delta functions. In some cases this can be done by careful 
counting of subspaces and dimensions [33, |35|, but it is not entirely trivial. In the 



FK representation, by contrast, everything is fully automated, but it is not sufficient 
to keep track of the connectivities of the sites in the top layer alone, as this layer will 
eventually need to be joined up to the bottom layer. Rather, it is necessary to keep 
track of the combined connectivities of the sites in the top and bottom layers: this 
method for handling periodic longitudinal b.c. will be sketched in Section [T4] and 



explained in detail in a subsequent paper [88 . 



3.3 Dimension of the transfer matrix 

Let us consider a strip of width m with free longitudinal boundary conditions. 
We want to know the dimension of the corresponding FK-representation transfer 
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matrix for different lattices (square and triangular) and different transverse boundary 
conditions (free and periodic). Each of these cases corresponds to a different class of 
allowed partitions V. The more stringent the restrictions on the allowed partitions, 
the smaller the dimension of the transfer matrix. 

We shall in general follow the notation of Stanley's 2-volume work Enumerative 
Combinatorics [ 117 , 118 1, where more detail on various integer sequences can be 
found. Another useful reference is the On-Line Encyclopedia of Integer Sequences 



119 1 . 



1) Let P m denote the set of all partitions of {1, . . . , m}. Its cardinality is given 
by the Bell numbers (or exponential numbers) B(m), which can be computed ||120 
Sections 6.1 and 6.3] from the exponential generating function 



E B (x) = y^Bim)^- = exp^-l) (3.42) 
/ — ' ml 

or from the remarkable formula 

B{m) = e" 1 ^^- (3.43) 

k=0 

(it is far from obvious at first sight that the right-hand side defines an integer!). For 
example, for m = 4 we have -8(4) = 15, corresponding to the partitions^ 

P 4 = {1, S±2, 5i3, #14, $23, #24> <^34> ^123j $124, ^134j <^234> ^1234> 

^12^34, S 13 S 2 4, ^Mfe} • (3.44) 

The Bell numbers have the asymptotic behavior 1|120| , Sections 6.1-6.3] [|121|, Section 
5.8] 1122], Examples 5.4, 5.10, 12.5 and 12.6] 

\ogB(m) = mW(m+l) - m + ™ +1 - - Hog[W(m+l) + 1] - 2 + 0(^H) 

W (m + 1) m 

(3.45a) 

TYl 

= mW(m) - m + ~ \ \og\W(m) + 1] - 1 + o(l) (3.45b) 

m log log m m /m(loglogm'' 2 
= m log m — m log log m — m + — + + O 



logm logm \ (logm) 2 

(3.45c 

as m — > oo, where W(z) denotes the unique real solution of we w = z (this is the 
Lambert W function [|123j| ) and has the (very slowly) convergent expansion in powers 

14 Henceforth we shall abbreviate delta functions 5(ai,(j 3 ) by (S 13 , etc. Moreover, we shall usually 
abbreviate partitions V by writing instead the corresponding product v-p of delta functions [cf. 
(^27))]: e.g. in place of V = { {1,3}, {2, 4}, {5} } we shall write simply V = S 13 5 24 ,- 
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of 1/ log z and log log z / log z: 

W(z) = hgz - log log z + V V -l n+1 n J 3.46 

/c! (log*)" 

where the s(n, £) are the Stirling numbers of the first kind. [The expansions ( |3.45| a,b) 
thus yield much better numerical approximations than ( |3.45) c) .] The Bell numbers 
for 1 < m < 14 are shown in Table 

The Bell numbers can also be written as 

B(m) = J2 S ( m ' k ) ' ( 3 - 47 ) 

k=l 

where the Stirling number of the second kind S(m, k) is the number of ways of par- 
titioning the set {l,...,m} into k nonempty subsets (i.e., the number of ways of 
placing m labeled balls into k indistinguishable boxes, with each box containing at 
least one ball). 

2) If the graph G is planar, then only non-crossing partitions can occur. The 
number of non-crossing partitions of {1, ... , m} is given by the Catalan number^ 

Cm = = -JL^ ( . (3.48) 



m!(m + l)! m+1 \ m 
Thus, for m = 4 we have C m = 14, corresponding to the partitions 

Pnc,4 = {1, <5l2, <5i3, 5\4, 623, #24, ^34, ^123, ^124, ^134, ^234, ^1234, 

^34,^23}, (3.49) 

since the only crossing partition of {1,2,3,4} is 613624. The Catalan numbers have 
the asymptotic behavior 

C m = 4 m m- 3 / 2 7r- 1 / 2 [l + 0(l/m)] (3.50) 

as m — > 00. For 1 < m < 14 they are shown in Table ||. 

3) For the zero-temperature antiferromagnetic model (v = —1), we can work in 
the basis w-p defined by ( |3.29| ), so that partitions with nearest neighbors in the same 
block are also forbidden. Let us consider first the case of free transverse boundary 
conditions. The number of non-crossing non- nearest-neighbor partitions of {1, ... , m} 
is given by the Motzkin number M m „i [|124|, Proposition 3.6] |125|1 .P| Here M n is the 



15 Stanley [118, Exercises 6.19 and 6.25] gives 66 com binatorial interpretations and 9 algebraic 
interpretations of the Catalan numbers C n . See also [118, Exercise 6.24]. 

16 Warning: Several references use the notation m n to denote what we call M n ; and one reference 
[ 126 1 writes M n to denote a different sequence. 
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number of ways of drawing any number of nonintersecting chords among n points on 
a circle^]; it is given by the explicit formula 

\n/2\ 

M « = E [Ij ) c > (3 - 51) 

and satisfies both the linear recursion relation 



2n + 1 3n — 3 

M n = — -M n _i + -M n _ 2 + 5 n0 (3.52) 

n + 2 n + 2 

and the nonlinear recursion relation 

71-2 

M n = Mn., + J2 M kM n -2-k + SnO (3.53) 

with initial condition M n = for n < (see e.g. ||127| , |126| , |128| , |129|| ). The generating 
function M(x) = X]^Lo M n x n is 



M(x) = ^ . (3.54) 



For example, for m = 4 we have M 3 = 4, corresponding to the partitions 

r nc+nnn,4 = {1,^13,^14,^24} • (3-55) 

The Motzkin numbers have the asymptotic behavior 

Mn = 3 n n~ 3 / 2 p^[l + 0(l/n)] (3.56) 

2^71 

as n —>■ oo. The transfer matrix for a triangular-lattice strip of width m with free 
boundary conditions at zero temperature has dimension 

TriFree(m) = M m _x . (3.57) 

The numbers TriFree(m) = M m _i for 1 < m < 14 are shown in Table 0. 

4) Let us now consider the zero-temperature antiferromagnetic model with pe- 
riodic transverse boundary conditions. The number of non-crossing non-nearest- 
neighbor partitions of {1, ... , m} when it is considered periodically (i.e. when 1 and 



17 This interpretation of M n arises in enumerating the "non-magnetic connectivities" of a loop 
model on the square lattice 114, pp. 1436-1437]. See Stanley [118, Exercises 6.38 and 6.46(b)] for 
14 combinatorial interpretations of the Motzkin numbers M n . 
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m also are considered to be nearest neighbors) is given by [129, Section 3.2, family 
R2]0 

' fo] "'" ' (3.58) 



d„ 



Rr, 



for m > 2 



where the Riordan numbers (or Motzkin alternating sums) R m [ 127] , |126| , |129|| P^| are 
defined by R = 1, Ri = and 



TO— 1 



E(-d 



m-fc-l M 



for m > 2 



k=0 



and satisfy the linear recursion relations 

Rm = —Rm-l + ^m-l + ^mO 

2m — 2 3m, 

R>m. 



m + 1 

and the nonlinear recursion relation 

m— 1 



2m — 2 3m — 3 

Rm-l + —rRm-2 + »m0 



m + 1 



(3.59) 

(3.60) 
(3.61) 



Rr, 



TO— 1 — 



+ 



fc=0 



with initial condition i? m = for m < (see e.g. ||127| , |126| , |129 
function R(x) = J2m=o Rm% m is 



R{x) 



X 



(3.62) 
The generating 

(3.63) 



2x(l + x) 

For example, for m = 4 we have d^ = R^ = 3, corresponding to the partitions 

* nc+nnnCyl,4 = {l,5i3,5 24 } (3.64) 
The Riordan numbers have the asymptotic behavior 



d. 



Rr, 



3 TO m -3/ 2 3v2 [1 + 0(1/m)] 



7T 



(3.65) 



as m 



oo. The numbers d m for 1 < m < 14 are shown in Table 



18 Let d m be the number of non-crossing non- nearest-neighbor partitions of {1, ... , m} when it is 
considered periodically. We have di = di = 1 and M m -\ — d m -i + d m for m > 3. [Proof: For 
to > 3, a non-crossing non-nearest-neighbor partition of {1, . . . ,m} (considered linearly) either has 
1 and to in the same block or it doesn't. There are d m partitions of the latter type. And in any 
partition of the former type, we can consider 1 and to to be amalgamated into a single site 1' that 
is a neighbor of both 2 and m — 1; so there are d m _i partitions of this type.] The formula (3.59) for 
d m then follows by induction. 

19 In most of the literature (e.g. 127, |126[ ) these numbers are called 7 m . We have adopted the 
recent proposal of Bernhart 1 129 to name them after Riordan |127|] and denote them R m . 
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5) Further restrictions arise from spatial symmetries. For instance, if we consider 
the square-lattice strip with free transverse boundary conditions, then reflection with 
respect to the center of the strip is a symmetry of the system. We can therefore 
define equivalence classes of non-crossing non-nearest-neighbor partitions modulo re- 
flection.0 The dimension SqFree(m) of the transfer matrix is then given by the 
number of these equivalence classes. Since each equivalence class contains either one 
or two partitions (and some contain only one), we clearly have 

— - — — - — - < SqFree(m) < TriFree(m) . (3.66) 
For m = 4, for example, we have three equivalence classes: 

PsqFree,4 = {1, &L3 + ^24, 8u} ■ (3.67) 

We have been unable to compute an explicit formula for SqFree(m), or to find any 
known integer sequence that corresponds to SqFree(m). The asymptotic behavior is 
clearly the same as that of TriFree(m) within a factor 2, hence of order 3 m m~ 3 / 2 . 
(We conjecture that SqFree(m) ~ TriFree(m)/2, since "most" partitions are asym- 
metric. The data in Table |2| strongly suggest that this conjecture is true: the ratios 
SqFree(m)/TriFree(m) are roughly decreasing in m, albeit with some even-odd os- 
cillation, and seem clearly to be approaching 1/2. Indeed, at m — 14 the ratio has 
already reached ~ 0.5032.) 

6) For the square lattice with periodic transverse boundary conditions, both re- 
flections and translations are symmetries. We therefore define equivalence classes of 
non-crossing non-nearest-neighbor partitions modulo reflections and translations and 
the corresponding number SqCyl(m) of equivalence classes.^ Since each equivalence 
class contains at most 2m partitions (and some contain less), we clearly have 

^ < SqCyl(m) < d m . (3.68) 
2m 

For m = 4, there are only two such classes: 

PsqCyl.4 = {Ml3 + <M • (3.69) 

We have been unable to compute an explicit formula for SqCyl(m), or to find any 
known integer sequence that corresponds to SqCyl(m). The asymptotic behavior is 
clearly the same as that of d m within a factor 2m, hence the leading behavior is 
~ 3 m . (We conjecture that SqCyl(m) « d m /(2m), since "most" partitions are asym- 
metric. The data in Table |2| strongly suggest that this conjecture is true: the ratios 
2mSqCyl(m)/<i m are roughly decreasing in m, albeit with some even-odd oscillation, 

One could also drop the non-nearest-neighbor condition; this would be relevant at nonzero 
temperature. 

21 One could also drop the non-nearest-neighbor condition; this would be relevant at nonzero 
temperature. 
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and seem to be approaching a value near 1. At m = 14 the ratio has already reached 
« 1.14.) 

7) For the triangular lattice, reflection is not a symmetry; but if we have periodic 
transverse boundary conditions, then translations are symmetries. We therefore de- 
fine equivalence classes of non-crossing non-nearest-neighbor partitions modulo trans- 
lations, and the corresponding number TriCyl(m) of equivalence classes.^ Since each 
equivalence class contains at most m partitions (and some contain less), we clearly 
have 

— < TriCyl(m) < d m . (3.70) 
m 

For m — 4, there are again only two such classes: 

P T nCyl,4 = {Ml3 + 5 24 }- (3-71) 

Indeed, we have found that TriCyl(m) = SqCyl(m) for m < 7; only for m > 8 does re- 
flection symmetry impose additional constraints beyond those imposed by translation 
symmetry, so that TriCyl(m) > SqCyl(m). We have been unable to compute an ex- 
plicit formula for TriCyl(m), or to find any known integer sequence that corresponds 
to TriCyl(m). The asymptotic behavior is clearly the same as that of d m within a fac- 
tor m, hence the leading behavior is ~ 3 m . (We conjecture that TriCyl(m) « d m /m, 
since "most" partitions are asymmetric. The data in Table ^ strongly suggest that 
this conjecture is true: the ratios mTriCyl (m)/d m are roughly decreasing in m, albeit 
with strong even-odd oscillation, and seem clearly to be approaching 1. At m — 14 
the ratio has already reached ~ 1.01.) 

Although the dimension of the transfer matrix is given by SqCyl(m), SqFree(m), 
TriCyl(m) or TriFree(m), the basis with respect to which the vectors are expressed 
is considerably larger, namely the set P n c,m of all non-crossing partitions (or P m 
in case G is non-planar). Thus, the vectors produced at intermediate stages of the 
computation can be rather long; this is the main limiting factor in our numerical 
work. In order to save memory, we use the following tricks: 

• Label the partitions by an integer variable k — 1, . . . , B(m). (For large widths 
m, however, it is vastly more efficient to include only the smaller set of the 
non-crossing partitions k — 1, . . . , C m |I6[ |115|| ; this approach will be taken in 
future work in collaboration with Jesper-Lykke Jacobsen |SB|, B7J-) 

• Represent the vectors in sparse-vector format (i.e. representing explicitly only 
nonzero coefficients). This is automatic in Mathematica. (If, however, we 
consider only the set of non-crossing partitions, then typical intermediate vec- 
tors are dense, so that sparse- vector format is actually a disadvantage.) 

22 One could also drop the non-nearest-neighbor condition; this would be relevant at nonzero 
temperature. 
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The coefficients of the chromatic polynomial are extremely large integers, which 
grow rapidly with the width and length of the strip. We can reduce the magni- 
tude of these coefficients (and therefore save memory) by performing a change 
of variables u = q — qo, where go nes & t or near the barycenter of the roots. 
This barycenter lies at |i£|/|V|, where \E\ (resp. |V|) is the number of edges 
(resp. vertices) in the graph G; that is, it lies at half the average coordination 
number. It is therefore convenient to take qo = 2 (resp. go = 3) for the square 
(resp. triangular) lattice; this vastly reduces the size of the coefficients of the 



chromatic polynomial, as was noted already by Baxter [55 



Remark. As explained in Section |3.1| , when computing the transfer matrix for a 
triangular-lattice strip of width m and cylindrical boundary conditions, we use the 
following technical trick [53]: To take account of the diagonal bond joining the sites 
m and 1, it is convenient to consider a triangular strip of width m + 1 and, at the 
end of the computation, identify the spins at sites m + 1 and 1. This means that the 
number of partitions arising in intermediate steps of the computation is not C m , but 
C m+ \. However, the dimension of the transfer matrix is still TriCyl(ra). 



4 Eigenvalue Crossing and Isolated Limiting Points 

4.1 Computation of eigenvalue-crossing curves 

The limiting curve B is the locus of points in the g-plane where there are two or 
more dominant eigenvalues. Our approach is to compute first the locus of q values 
where there are two or more equimodular eigenvalues, dominant or not; we then check 
the corresponding eigenvalues one-by-one for dominance. We have used two methods 
to compute the locus of equimodularity: the resultant method, and the direct-search 
method. 



4.1.1 The resultant method 

The resultant of two polynomials P{x) = A Y[fLi( x ~ x i) an d Q{y) = B Y\^ = i{y— yi) 
is defined to be the product of all the differences of roots, scaled suitably by the two 
leading coefficients: 

M N 

Res(P, Q) = A N B M J] JJ(x, - y 3 ) . (4.1) 

i=l j=l 

Thus, the resultant vanishes if and only if P and Q have at least one root in common 
(or one or both of the leading coefficients vanishes). It is a nontrivial fact that the 
resultant Res(P, Q) can be expressed as an (M + N) x (M + N) determinant involving 
the coefficients of P and Q; it is not necessary to know explicitly the roots {x^ and 

{.'/,}• 
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Likewise, the discriminant of a polynomial P(x) = A Y\ i= defined as 

Disc(P) = A^Hixi-xj) 2 = (-l^-^A^Hixi-xj) . (4.2) 

i<j i^j 

(For a quadratic P(x) = ax 2 + bx + c, we have Disc(P) = b 2 — 4ac, which agrees with 
the definition used in high-school algebra.) It is not difficult to show that 

Disc(P) = (-l) M(M - 1)/2 A- 1 Res(P,P') , (4.3) 

so that the discriminant vanishes if and only if P has at least one multiple root.0 

For further information on resultants and discriminants and algorithms for com- 
puting them, see e.g. ||131| , Chapter 3]. 

Consider now the characteristic polynomial of the transfer matrix T(q): 

dimT 

P(X,q) = det[A/ - T(g)] = J] ( A " A ^)] > ( 44 ) 



i=l 



where {Aj(g)} are the eigenvalues of T(q). Consider next the polynomial 

P B (X,q) = P{e iQ \q). (4.5) 

P and Pq are polynomials in A whose coefficients are polynomials in q (and in e ), and 
they have a root in common if and only if T(q) has eigenvalues Ai and A2 satisfying 
Ai = e i<9 A 2 . (Note that, in addition to the desired case of an equimodular pair of 
eigenvalues, this also occurs whenever T(q) has a zero eigenvalue.) We can therefore 
compute the locus of equimodularity by sweeping over a closely spaced set of points 
9 G (0, 7r], and computing for each 9 the roots of 

R e {q) = Res A (P,P e ) , (4.6) 

which is a polynomial in q and e l0 .f^\ 

Of course, 9 = is a special case: here we are looking for the multiple roots of P, 
which can be located by computing the zeros of the discriminant 

*< S ) = Res> ( p,n = (4 . 7) 

where P'(X,q) = dP(X,q)/dX. As we shall see, the zeros corresponding to 9 = are 
very important because they correspond to the endpoints of the curves of equimod- 
ularity. 



23 Lang [130, p. 204] warns that many books (including his own first edition!) contain sign errors 
(i.e. fail to include the factor (— 1) M ( M_1 )/ 2 ) either in the definition of the discriminant or in the 
formula relating it to the resultant. 

24 This method was presumably known already to Beraha, Kahane and Weiss: see |58| footnote 

21. 
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After finding a set of q values for which T(q) has a pair of equimodular eigenvalues, 
we check them one-by-one by solving the characteristic equation P(\, q) = and 
testing whether the pair of equimodular eigenvalues is dominant or subdominant. 
The limiting curve B of partition-function zeros corresponds only to the crossing in 
modulus of dominant eigenvalues; nevertheless, it is sometimes of interest to depict 
the loci of subdominant eigenvalue-crossing as well. Whenever we do so in this paper, 
we shall draw the dominant eigenvalue-crossing curves using solid black lines and the 
subdominant curves with dashed red lines. 

In practice, we first use Mathematica to compute Re(q) symbolically as a poly- 
nomial in q and e ld with integer coefficients. To minimize the effect of round-off errors 
in the subsequent computation, we use instead of 9 the real variable t = tan(#/2) e 
(0, oo] defined by 

,-fl 1 + it , . 

e = — t • < 4 - 8 > 

and we always choose t to be a rational number. In this way, e td is a complex rational 
number whose modulus is always exactly equal to 1, and the polynomial Re{q) has 
complex rational coefficients. This allows us to take advantage of arbitrary-precision 
polynomial root-finders such as Mathematical NSolve or (better yet) the package 
MPSolve 2.0 developed by Bini and Fiorentino |32|, |133j.pl 



The drawback of the resultant method is that the degree of the resultant poly- 
nomial Re(q) grows very rapidly with the width of the lattice strip; moreover, the 
coefficients in Rg(q) also grow very rapidly (even if we use the variable u = q — qo). 
This limits in practice the widths we can study: L x < 6 for free transverse boundary 
conditions and L x < 8 for periodic transverse boundary conditions. 

Remark. In some cases (notably with periodic and twisted-periodic boundary condi- 
tions in the longitudinal direction |26], ^7[ |29|, [H], [32|, |33|), the amplitude correspond- 
ing to one or more eigenvalues can vanish identically. If this occurs for one of the 
dominant eigenvalues, blind application of the foregoing procedure can lead to erro- 
neous results. As a safeguard against this phenomenon, we numerically computed the 
corresponding amplitudes for a selected value of q, by numerically diagonalizing the 
transfer matrix and rotating the corresponding vectors on left and right. It suffices 
to find at least one value of q for which none of the amplitudes vanishes. In all the 
cases considered in this paper, an identically-vanishing amplitude never arises. 

Example. In the simple special case of a 2 x 2 transfer matrix 

r(9) - ( t Z ) (4 - 9 » 



25 MPSolve 2.0 is much faster than Mathematical NSolve for high-degree polynomials (this 



is reported in |133 |, an d we confirm it); it gives guaranteed err or bo unds for the roots, based on 



available [ 132 



rigorous t heor ems 1133] ; its algorithms are publicly documented [133]; and its source code is freely 



27 



with left and right vectors 



u ( f } q ] ) (4.10a) 



9(q) 

1 

o 



(4.10b) 



and partition function 

Z n = u T T n ~ l v (4.11) 

for a strip of length n (where T denotes transpose), we can obtain explicit expressions. 
We have 

Z n (q) = a_(g)A_(g) B - 1 + a + (g)A + (g)^ 1 , (4.12) 
where the eigenvalues A± and their corresponding amplitudes a± are 

A±(g) = \ (taT(g) ± ^t^T{q) - 4detT(g)) (4.13a) 

= ^(A(9)±V^)) (4.13b) 

and 

■ ^ |/fe,± J&)- (414b) 

Note that even though the eigenvalues and the amplitudes are non-polynomial alge- 
braic functions of q, the partition function Z n (q) is always a polynomial in q. From 
these equations it is not difficult to prove the recurrence relation 

Z n+2 {q) - [txT{q)]Z n+l {q) + [det T(q)]Z n (q) = 0. (4.15) 

According to the Beraha-Kahane- Weiss theorem (Theorem |2.1| ), the limit points 
as n — > oo of the partition-function zeros are of two types: the isolated limiting 
points, which occur where one eigenvalue is dominant and its amplitude vanishes; 
and curves of non-isolated limiting points, where there is a crossing in modulus of 
two or more dominant eigenvalues. In our simple case, the limiting curves are given by 
the condition |A_(g)| = |A + (g)|, which means that P 2 (q)/Pi(q) 2 should be a negative 
real number (call it — t 2 ), or in other words 

tr 2 T(g) = 4?detT(g) with < t = < 1 . (4.16) 

1 — (~ t 



An identical result is obtained using the resultant method ( |4.6|) /(477): we have 

P(\,q) = A 2 - [trT(g)]A + detT(g) (4.17) 
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and hence 



R e (q) = (l-e i0 ) 2 [detT(q)} [(1 + e ie ) 2 detT(q) - e ie tr 2 T(q)] (4.18a) 

R (q) = 4detT(g) -tr 2 T(g) = -P 2 {q) (4.18b) 

We see that (f4. 18a ) vanishes precisely on the curve ( 4.16| ): it suffices to insert (|4.8j ) 
and note that 

(l + e ie ) 2 4 

— P- = T^—, = 4t (4.19) 
e* e 1 + t 2 K J 

[equivalent ly, t = cos 2 (6 l /2)]. 

The isolated limiting points of zeros are given by the condition that the ampli- 
tude of the leading eigenvalue vanishes. We first compute the product of the two 
amplitudes: 

MffK(?) = ( 4 - 2 o) 

where 

PM = f{q)2P2{q l - PM2 = zmzM - zm 2 ■ (4.2i) 

[Note that Q4~20D / (|OID is just Lemma specialized to the case M = 2.] We now 
observe that a root of P4 corresponds to the vanishing of a_ (resp. a + ) in case f \fP~2 = 
P3 (resp. fy/Pz = —P3)', and A_ (resp. A + ) is dominant in case Re(Pi/v^) < 
(resp. > 0). It follows that a root of P 4 corresponds to the vanishing of the amplitude 
corresponding to the dominant eigenvalue in case 

Re^iM < Q (4 22) 

there. If Ps{q) = at a root of P4, then both amplitudes vanish there. 



4.1.2 Direct-search method 

When we were unable to apply the resultant method (i.e. for large lattice widths), 
or when we wanted to study in detail a small region in the g-plane (for any lattice 
width), we used a direct-search method to obtain the curves of equimodularity. The 
idea is to define a function that measures the difference between the moduli of the 
two dominant eigenvalues of the transfer matrix. Since an explicit expression of the 
eigenvalues as a function of q is usually not available, the eigenvalues are obtained 
numerically for each value of q. Then we extract the two eigenvalues of largest 
modulus and compute 

F(q) = |Ai(g)|-|A a (g)| (4.23) 

where u>Xi(q) ^ u\2{q) in lexicographic order (here u is some arbitrarily chosen 
nonzero complex number). Clearly F(q) vanishes if and only if the two dominant 
eigenvalues are equimodular; moreover, because generically the two dominant eigen- 
values will not satisfy Re[a;Ai(g)] = Re[cl;A2(q , )] precisely where they also satisfy 
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|Ai(g)| = ^2(9) |, they will generically not interchange roles at the eigenvalue crossing, 
and hence F(q) will be a smooth function of q there.0 We then search for the zeros 
of F(q) using a Newton method in the complex g-plane. 

Once a good approximation for a zero is found, we also compute the phase e ld = 
A1/A2 [or equivalently, the number t defined by (|4.8|)1. Knowledge of 6 is very useful 
in trying to understand the topology of the limiting curve. 



4.2 Qualitative structure of eigenvalue-crossing curves 

In this section we want to discuss the types of qualitative behaviors that can 
arise when studying the eigenvalue-crossing curves. Recall first [|134| , Chapter 2] that 
the eigenvalues Aj(g) of the transfer matrix T(q) are analytic functions of q except 
possibly where two or more eigenvalues collide. Indeed, this can be seen by expanding 
the characteristic polynomial P(X, q) = det[AJ — T(q)] around a root (A , go): 

P(X,q) = a(X-X ) + b(q-q ) + c(X-X ) 2 + d(q-q ) 2 + e(X-X )(q-q ) + ... (4.24) 

Provided that A is not a multiple eigenvalue of T(g ), we have a = (dP/dX)(X , q ) 7^ 
0; the implicit function theorem then guarantees that in a neighborhood of g = go 
there exists an analytic function X(q) solving P(A(g), q) — with A(go) = Ao, and it 
has the convergent expansion 

A(g) = Ao - -(g - g ) + — " a J ~ — (g - g ) 2 + 0((g - go) 3 ) ■ (4.25) 

a a 6 

If, on the other hand, A is a fc-fold eigenvalue of T(g ) with k > 2, then A(g) can 
have an l th -Toot branch point at go for any I < k. More precisely, near g = g the 
eigenvalues divide into groups of cardinalities h, ■ ■ ■ ,1m (with h + ■ • ■ + Im — dimT) 
such that the eigenvalues of the i th group are the U distinct determinations of an 
analytic function of (g — go) 1 ^', i-e. they have an Z^-root branch point at go ||134, 



pp. 65-66]. For example, for k = 2 we have a = and c 7^ 0: if b 7^ 0, then a pair of 
eigenvalues A±(g) bifurcate off A with a square-root branch point, 



A±(g) = Ao ± (--X 2 (g - g ) 1/2 - ^(g - go) + 0((g - go) 3/2 ) ; (4.26) 
while if b = 0, then A±(g) are analytic functions in a neighborhood of g , 

A ± (g) = A + ~ 6±V 2 e c 2 ~ 4cC? (g - g ) + 0((q - g ) 2 ) . (4.27) 



26 There is nothing special about lexicographic order; virtually any order will do, except ordering 
by moduli. In practice, we used Mathematical default ordering, which is slightly different from 
lexicographic. Note also that any nonzero u> £ C will do; but when studying an eigenvalue crossing 
on the real q axis (for which Ai and A2 will be complex conjugates), it is advantageous to choose w 
non-real. 
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4.2.1 Crossing of two simple eigenvalues 



Generically, the equimodularity of two eigenvalues defines an analytic curve in the 
complex g-plane, along which the parameter 9 (or t) varies smoothly. To see this, 
suppose that at q = q we have a pair of simple (i.e. non-multiple) eigenvalues Ai,o 
and A2,o that have equal modulus (| Ax,o| — I ^2,0 1 7^ 0); they satisfy Ai,o = e l8 X 2 ^ with 
6 7^ (mod 27r). Each of these eigenvalues then extends to a single- valued analytic 
function of q in a neighborhood of q — qo, as in ( 4.25) : 



Aj,o [q 

Cti 



qo) + 



dihei - afdi - bfc 
a? 



iq-qof + O((q-q ) s ) (4.28) 



for i = 1,2. Their ratio is then 
AiO?) 



jo 



A 2 (g) 



e ie [l + pe i %-q ) + O({q-q 



(q - q ) + 0{{q - q ) 



(4.29a) 
(4.29b) 



Provided that p 7^ 0, the equimodularity locus | Ai (g) / A2 (g) | = 1 defines near q = q 
an analytic curve that passes through q = qo at angle —0 ± tt/2. 

If, however, p = [i.e. the term in ( |4.29aj ) that is linear in q — q vanishes], 
then the equimodularity locus can have multiple points. (A sufficient though not 
necessary condition for this to occur is for the linear terms to vanish in each eigenvalue 
separately, i.e. b\ = 62 = 0.) If the first nonvanishing term in ( 4.29a|) is the one of 
order (q — qo) k , then we have a fc-fold multiple point in the sense of algebraic geometry 

135], Section 6.2]: that is, the equimodularity locus \Xi(q) / X 2 (q)\ 

i = qo an analytic image of the set Re z k = near z 



1 defines near 
0. In particular, this set can 



be interpreted locally as k analytic curves crossing at angles ir/k. 

Remark. Let us consider the partition function Z n (q) near a point go where there 
are exactly two dominant simple eigenvalues: 



Z n (q) = a 1 (q)X 1 (q) n + a 2 {q) X 2 {q) n 



(4.30) 



where the dots indicate the contributions of subdominant eigenvalues. Let us assume 
for simplicity that Oi(g ), «2(<?o) 7^ and insert the expansions 



log 
log 



Ai(g) 

Hq) 



a 2 (q) 



2vr# + A(q-qo) + 0((q - qof) 



2mB + C{q-qo) + 0((q - q ) 2 ) 



where ip = 9/2n E R and A,B,C E C. [From (4.29) we have A 
(bi/aiXifi).] Then ( |4.30| ) can be written as 



pe 



i(f> 



(4.31) 

(4.32) 
(& 2 /ct2A 2 ,o) - 



Z n {q) = 2t[a 1 (q)a 2 (q)X 1 (qrX 2 (qr] 1/2 x 

An + C 



sinh 



mvmj} + B) + 



(q-q ) + O((q-q ) 



(4.33) 
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which agrees with |55|, eqn. (3.3)] by trivial renaming of variables. When n is large, 
( |4.33| ) has zeros at 



q = qo - -r[k + B + R{nip)) + 0{l/n 2 ) (4.34) 

where k G Z is an arbitrary integer and R(x) = x — \_x\ is x modulo 1. Therefore, we 
see that near a regular point go of the limiting curve B, the finite-volume partition 
function Z n (q) has a sequence of evenly spaced zeros [spacing = 27r/(|A|n)] lying 
parallel to B. These zeros lie a distance 27r|Im5|/(|A|n) away from B, hence converge 
to it generically at rate 1/n. 

In addition to this generic behavior, there are other features exhibited by the 
limiting curves that are worth studying in detail: 



4.2.2 Endpoints (collision of two eigenvalues) 

Suppose that A ^ is a two-fold eigenvalue of T(q ). Then generically (i.e. if 
b 7^ 0) a pair of eigenvalues \±(q) bifurcates off A with a square-root branch point, 
as in ( |4.26D . The ratio of the two eigenvalues is 



= 1 + A(q - q ) 1/2 + y (? - go) + 0((q - q G f 2 ) (4.35) 

where A = (2/A )(-6/c) 1/2 ^ 0. Setting 

= e* = ^ (4.36) 
A_(g) I -it ( J 

with t real, we see that the equimodularity locus is given by 

q = qo ~ ^t 2 + 0{t') . (4.37) 

Only even powers of t appear in this expansion, since A + and A_ interchange roles as 
we go around the branch point, and A + /A_ = e l6 if and only if A_/A + = e~ ie . Writing 
A = pe 1 ^, we see that the equimodularity locus is an analytic curve ending at q = q 
and tangent there to the ray arg(g — q ) = n — 2(f>. 

Because the two eigenvalues are equal at q = q , the parameter 6 (or t) takes the 
value at endpoints. Each such endpoint corresponds to a simple root of the t = 
resultant Ro(q)'- indeed, ( |4.26| ) implies that for q near q we have A + — A_ ~ (q— qo) 1 / 2 
and hence Ro(q) ~ (A + — A_) 2 ~ q — q . For this reason, the simple zeros of the t = 
resultant play an essential role when we try to determine with high accuracy the 
topology of the limiting curve. 

By contrast, in the non- generic case 6 = 0, where the eigenvalues A±(g) are an- 
alytic functions in a neighborhood of go [cf. ( [4.271 )1, the formulae (|4.28|) /(4.29) of 



the preceding subsection apply with 6 = 0. In this case, the root q = go is a 
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double zero of the t = resultant Ro(q), as A + — A_ ~ q — go by ( 4.271) , so that 
Ro(q) ~ (A+ - A_) 2 ~ (g - g ) 2 . 

Example. Consider a two-dimensional transfer matrix as in ( [4.91 ). Collision of the 
two eigenvalues ( [4.13|) occurs whenever p2{q) = 0, i.e. whenever tr 2 T(g) = 4detT(g). 
Suppose this occurs at q = q , with A = A±(g ) = |trT(g ), and let us expand the 
eigenvalues (|4.13|) around q = g : 



A± = ^trT(go) ± Ci(g-g ) 1/2 + C 2 (q - q ) + 0((g - g ) 3/2 ) 



where 



C 2 



1 dtiT(q) 
2 



(idetT(g) 



9=90 



(/(/ 



1/2 



9=90 , 



(if/ 



9=90 



(4.38) 

(4.39a) 
(4.39b) 

(4.40) 



The ratio of the two eigenvalues is therefore 

^±|| = 1 + A{q-q Q f' 2 + ^(g-fc) + O((q-q f 2 ) 

with A = 4d/[trT(g )]. 

Remark. Let us return now to the general case, and consider the partition 
function Z n (q) close to an endpoint q = qo where two dominant eigenvalues collide: 

Z n {q) = a + (q)\ + (q) n + a_(q)\_(q) n + ... (4.41) 

where the dots indicate the contributions of subdominant eigenvalues. From ( |4.26|) 
[cf. also HD/dH)] we have 

A+(g) 



log 
log 



A-((7) 
a+(q) 



A{q - g 



,V2 



a_(g) 



0(q - g ) 



2mB + C(g-g ) 1/2 + 0{q - g 



(4.42) 
(4.43) 



where A = 2(—b/c) 1 ^ 2 /Xo- Then ( |4.41|) can be written as 
Z n (q) = 2z [a + (q) a^q) X + (q) n X4q) n } 1/2 x 

An + C 



sinh 



inB + 



(g-g ) 1/2 + O(g-go) 



When n is large, the zeros of (|4.44| ) are given by 

4vr 2 (A; + B) 2 



go - 



A 2 n 2 



+ 0(l/n 3 



(4.44) 



(4.45) 



where k e Z, in agreement with ( |4.37|) . Therefore, the convergence close to an 
endpoint is generically of order 1/n 2 , in contrast with the 1/n convergence near a 
regular point of the limiting curve [cf. ( |4.34 )j. 
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4.2.3 Crossing of three simple eigenvalues (T points) 



Suppose that at g = go we have three dominant simple eigenvalues (call them 
-^1)^2, A3): then three smooth curves of equimodularity, corresponding to the three 
pairs of eigenvalues A1/2, A1/3, A2/3, pass through q = q . One half of each curve of 
equimodularity corresponds to a dominant pair of eigenvalues, while one half corre- 
sponds to a subdominant pair; therefore the locus of dominant equimodularity looks 
like a T, so we call these crossings T points (see Figure H). They are complex ana- 
logues of "triple points" in the thermodynamic sense (but they are not multiple points 
in the sense of algebraic geometry: see Remark 2 below). 

More precisely, the eigenvalues Xi(q) for % = 1,2,3 define single- valued analytic 
functions (|4.28|) in a neighborhood of go, so that their ratios are 



Ai(?) 


Ai,o 


1- 


( a x 


a.2 


A 2 (g) 


A2,0 


VMi,o 


02X2,0 


Ai(g) 


Ai,o 


1- 


( ax 


a 3 


A 3 (g) 


A3,0 


VMi,o 


^3A3,o 


A 2 (g) 


A 2 ,o 


1- 


( a 2 


a 3 


A 3 (g) 


A3,0 




^3A 3j o 



(g-g ) + O((g-g ) 2 ) 
(q - qo) + 0((q - g ) 2 ) 
(q - 9b) + 0((q - g ) 2 ) 



(4.46a) 
(4.46b) 
(4.46c) 



Provided that the coefficients of g — go in ( 4.46 a-c) are not collinear, there are six 
possible phases for g — g that render one of the ratios unimodular to leading order in 
q — qo', three of these correspond to a dominant crossing, while the three complemen- 
tary phases correspond to a subdominant crossing. If, on the curve where |Aj/Aj| = 1, 
we define 



Ai(g) 



(4.47) 



then the phase %(g) varies smoothly along this curve. No pair of these three phases 
obeys in general any relation, even as q — > go; but at the T point we obviously have 
the relation 

^12 + ^23 + ^31 = 2vm (4.48) 

for some integer n. If we take —it < 9ij < n, then n must be —1, or +1. In particular, 
the absolute values of the #y (which are what we actually calculate in practice) must 
satisfy either \6i 2 \ + 16*23 1 + I #31 1 = 2ir or some permutation of 1 6*12 1 + | ^23 1 — |#3i| = 0. 
In practice, one way of detecting a T point is by noticing the discontinuity in 9 as we 
pass from one piece of the limiting curve to another. 

Remarks. 1. We frequently find that the limiting curve B contains an arc starting 
at an endpoint (t = 0) and ending at a T point (t = to)-0 ^ we are a ^^ e ^° compute 
the t = resultant, then all endpoints can be detected, so that the corresponding 



27 Similar features were found in the complex-temperature zeros of the two-dimensional spin-s 
Ising model [ 136 1 , in the complex-temperature zeros of the Potts model on the hexagonal, Kagome 
and triangular lattices fl5|, [l6| , and in the complex-g zeros of the Potts antiferromagnet on some 
families of strip graphs 
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arcs can be detected as well. But using the direct-search method, a short arc can 
easily be overlooked: for since the parameter t grows smoothly as we move along the 
arc towards the T point, a short arc corresponds to a small value of t ; and if to 
is smaller than our step size, we may fail to detect the discontinuity in t along the 
other two arcs merging at the T point (unless we are very lucky!). This means that 
in the cases where we were unable to obtain the zeros of the t — resultant (namely, 
widths 7p and 8p), we may have missed some such "close pairs" of endpoints and T 
points; so the lists of endpoints and T points reported in Sections |5.6| and [5]7| may be 
incomplete, and the counts reported in Table [5] must be understood as lower bounds 
on the true value. 

2. Rocek, Shrock and Tsai |17], p. 528 top, item (ii)] state erroneously that a T 
point is a "multiple point" in the technical sense of algebraic geometry, i.e. a point 
where two or more branches of the same irreducible algebraic curve meet. In fact, 
a T point corresponds to the crossing in modulus of three usually unrelated analytic 
functions. Our aim here is not to quibble about definitions, but to emphasize a radical 
difference in behavior: at a fc-fold multiple point the branches always intersect at 
angles n/k, while at a T point they can (and in general do) intersect at arbitrary 
angles. (By contrast, the crossing (4.29) with p = is a true multiple point, as 
correctly observed by Rocek et al. |T7|, p. 528 top, item (i)].) 



4.3 Computation of isolated limiting points of zeros 

According to case (a) of the Beraha-Kahane-Weiss theorem, the isolated limit- 
ing points of partition-function zeros correspond to points where there is a unique 
dominant eigenvalue whose amplitude vanishes. We locate such points by a two-step 
process: First we determine, using Lemma ^2, the (finite) set of q values where at 
least one amplitude vanishes. Then we check these q values one by one, by diagonal- 
izing T(q), rotating the left and right vectors u(q) and v(q), and checking whether 
the vanishing amplitude(s) corresponds to the dominant eigenvalue. In the special 
case dimT = 2, we can test dominance using the analytic criterion ( |4.22| ). 



5 Numerical Results for the Square-Lattice Chro- 
matic Polynomial: 

Free Transverse Boundary Conditions 

We have computed the transfer matrix T(q) and the limiting curves B for square- 
lattice strips of widths 2 < L x < 8 with free boundary conditions in the transverse 
direction. We have checked the self-consistency of our finite-lattice results using the 
trivial identity 

Z{m-p x n-p) = Z{n-p x mp) (5.1) 

for all pairs 2 < m, n < 8. 
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5.1 L x — 2p 

In this case the transfer matrix is one- dimensional, and the result is trivial: 

Z(2 F x n F ) = q(q - l)(q 2 - 3q + 3)"- 1 . (5.2) 

Since there is only one eigenvalue, there is obviously no crossing, hence B = 0. 
However, there are zeros for all n at q — 0, 1 (trivially) and for all n > 2 at q = 
(3±V3i)/2. 



5.2 L x = 3f 

In this case the transfer matrix is two-dimensional. The allowed partitions are 
given by P = {1, (J13}. In this basis the transfer matrix is equal toQ 

T(3 F ) = (V-5^ + 10,-8 t-M + ty (5 . 3) 
and the partition function is equal to 

T 

Z(3 F xnpj^fg-l)^ 1 ) ■ T(3 F )- 1 .(J) (5.4) 

where T denotes transpose. 

We can rewrite the above expression for the partition function as in ( J4.12p . The 
polynomials P±, P2, P3 and P4 entering the definitions of the eigenvalues and ampli- 
tudes (PD/(CT)/(PD|) are given by 

Pi(q) = (g 2 -3g + 5)(g-2) (5.5a) 
P 2 (q) = ( g 2_ 5g + 7 )( g 4_ 5g 3 + llg 2_ 12g + 8 ^ ( 55 ^ 

P 3 ( q ) = g(g_l)( g 4_ 6g 3 + 14g 2_ 15g + 8 ^ 

PM = q 2 (q-l)\q-2) (5.5d) 

The limiting curve B (see Figure 0) consists of three disjoint arcs.0 One of them 
crosses the real axis at go — 2, and is invariant under complex conjugation; the other 
two lie in the first and fourth quadrants, respectively, and are complex conjugates of 
each other. There are six endpoints: 

q w 0.5865699800 ± 1.1400627519 i (5.6a) 
q w 1.9134300200 ± 1.0979688996 i (5.6b) 

q = 5± w 2.5 ± 0.8660254038 i (5.6c) 



28 As was found a quarter of a century ago by Biggs and Meredith |111, p. 11]. 

29 This curve is also depicted in jl7| Figure 3(a)]. 
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These endpoints are the six roots of the t — resultant Ro(q) = —P2(q) given by 
(5.5b). 

The zeros of the amplitudes can be found by solving the equation P/t(q) = 0. 
There are two trivial zeros q = 0, 1, where both amplitudes vanish simultaneously; 
both these zeros lie in the region where the eigenvalue A_ is dominant. The non-trivial 
zero q = B 4 = 2 is a zero of a_, but not of a + [a + (q = 2) = 2]. However, this point 
happens to lie on an dominant-eigenvalue-crossing curve [X+(q = 2) = — A_(g = 2) = 
1], so it actually belongs to case (b) of the Beraha-Kahane- Weiss theorem. Indeed, 
from Table |3] we see that the first non-trivial real zero of the partition function seems 
to converge to q = 2, but at the slow (roughly 1/n) rate characteristic of non-isolated 
limiting points rather than at the fast (exponential) rate characteristic of isolated 
limiting points. 



5.3 L x = 4p 

The transfer matrix is three-dimensional. In the basis P = {l,<5i3 + ^24,^14} it 
can be written as 

T(4 F ) = 

g 4 -7g 3 + 21g 2 -32g + 21 2(g 3 - 6g 2 + Uq - 12) g 3 - 7q 2 + 19q - 20 
q-2 q 2 - 4g + 5 3-g 



-1 -2(g-2) q 2 -5q + 7 



(5.7) 



and the partition function is equal to 

/(?-l) 2 \ T 

Z(A F xn F ) = q(q-l) 2{q - 1) • T^p)"" 1 • I I . (5.8) 





The limiting curve B (see Figure ^) has three connected components: again, one 
crosses the real axis and is self-conjugate, while the other two stay away from the real 
axis and are a pair of mutually conjugate arcs.f^] This time, however, the component 
that crosses the real axis is rather complicated: there is a pair of T points at q ~ 
2.327 ± 0.9113 i, and there is a double point on the real axis at q w 2.2649418565. 

There are ten endpoints: 

q w 0.3254743549 ± 1.1048503376 i (5.9a) 

q w 2.0555822564 ± 1.5703029256 i (5.9b) 

q w 2.2283590792 (5.9c) 

q « 2.2823594125 ± 1.5512247035 i (5.9d) 

q w 2.3014157308 (5.9e) 

q w 2.6674264726 ± 0.7845284722 i (5.9f) 



:!() 



This curve is also depicted in mn Figure 3(b)]. 
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The horizontal segment emerging from the double point ends at the pair of real 
endpoints (|5.9|c,e) . 

Using Lemma [2.2| we find that the points where at least one amplitude vanishes 
are given by the zeros of the polynomial 

det D(q) = 2q 3 (q - l) 3 (g - 2) 2 {q 2 - 3q + l)(2g 3 - 13g 2 + 27q - 17) 2 . (5.10) 

The values q = 0, 1 are trivial zeros where all three amplitudes vanish simultaneously. 
At q = 2 the amplitude corresponding to the leading eigenvalue vanishes (as does 
one of the two subdominant amplitudes). The other roots q — (3 ± y/E)/2 and 
q pa 1.170516, 2.664742 ± 0.401127 i all correspond to the vanishing of a subdominant 
amplitude. Therefore, according to case (a) of the Beraha-Kahane-Weiss theorem, 
the isolated limiting points for this strip are q — 0, 1, 2. 

From Table |3|, we see that the first non-trivial real zero converges rapidly to the 
Beraha number _B 4 = 2. In addition, there are further real zeros (whose number 
increases with n) that tend to the segment [2.2283590792 . . . , 2.3014157308 . . .] of the 
limiting curve and in the limit n —>■ oo become dense on that segment. 

It is curious that (|5.10|) vanishes also at the Beraha number i?s = (3 + Vb)/2 and 
its conjug ate B* = (3-V5)/2, even though both of these correspond to the vanishing 
of a subdominant amplitude. 

5.4 L x = 5f 

The transfer matrix is seven-dimensional, and can be expressed in the basis P = 
{1, Si 3 + 5 35 , 5 2 4, ^14 + ^25, ^15, ^135, ^15^24}- We refrain from giving here the full transfer 
matrix; instead, we refer the reader to the MATHEMATICA file transfer 1 .m included 
in the electronic version of this article at xxx.lanl.gov. 

The limiting curve B (see Figure |5|) has five connected components. One of them 
crosses the real axis at go ~ 2.4284379020 and is a self-conjugate arc; there is a pair of 
mutually conjugate arcs; and finally, there is a pair of mutually conjugate components 
exhibiting T points at q pa 2.423 ± 0.1067z and q pa 2.291 ± 1.5611 

There are 14 endpoints: 



q ^ 0.1708973690 ± 1.0464583589 i (5.11a) 

q pa 1.9065720451 ± 1.9339587717 i (5.11b) 

q pa 1.9748200483 ± 1.9395387106 i (5.11c) 

q pa 2.3024178902 ± 1.6190810539 i (5. lid) 

q pa 2.3990745384 ± 0.8206408701 i (5. lie) 

q pa 2.4983799650 ± 0.8199051472 i (5.11f) 

q pa 2.7692051339 ± 0.7143320949 i (5.11g) 



Please note that the endpoints q pa 1.9065720451±1. 9339587717 z and q pa 1.9748200483± 
1.9395387106 i define a pair of small gaps in the limiting curve. Likewise, the curves 
emerging from the endpoints q pa 2.3024178902±1.6190810539 % and q pa 2.3990745384± 
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0.8206408701 i do not cross, but define another pair of very small gaps (see Figure |6] 
for detail). 

By Lemma |2]2|, the points where at least one amplitude vanishes are given by the 
zeros of the polynomial 

det D(q) = q\q - l) 13 (g - 2) V - 3q + l) 5 (q - 3) 7 P(q) 2 (5.12) 

where P(q) is a polynomial of degree 37 with integer coefficients that we report in 
the file transfer l.m; as far as we can tell, P(q) cannot be factored further over 
the integers. The values q = 0, 1 are trivial zeros where all seven amplitudes vanish 
simultaneously. (Two of the amplitudes have in fact multiple roots at q = 1.) At q = 2 
the amplitude corresponding to the leading eigenvalue vanishes (as do the amplitudes 
of five of the six subdominant eigenvalues). In addition, one complex-conjugate pair of 
zeros of P(q), namely q ~ 2. 2866147868±1. 0116506019 i, corresponds to the vanishing 
of a dominant amplitude (shown with a x in Figure ||). This is the first time we find 
a nonreal isolated limiting zero. The remaining zeros of det-D(g) correspond to the 
vanishing of subdominant amplitudes. 

The first non-trivial real zero (see Table 0) converges quickly to the Beraha number 
Bi = 2. The next real zero appears to be converging slowly (at a roughly 1/n rate) 
to the value go ~ 2.4284379020 where the limiting curve B intersects the real axis, 
which is smaller than the next Beraha number B 5 . The convergence to the complex 
isolated limiting points at q ~ 2.2866147868 ± 1.0116506019 i is quite again rapid. 
Indeed, if we select for each strip length n the zero closest to the limiting point, and 
fit the distance from the limiting point to an inverse power of n, the effective power 
seems to grow with n; this is compatible with the expected exponential convergence. 

It is curious that ( |5.12| ) vanishes also at the Beraha numbers B 5 = (3 + v^5)/2 
and B 6 = 3 — and hence also at the conjugate B£ = (3 — y/E)/2 — even though all 
of these correspond to the vanishing of a subdominant amplitude. 

5.5 L x = 6p 

The transfer matrix is 13-dimensional; it can be found in the Mathematica file 
transfer 1 .m. 

The limiting curve B (see Figure [7]) has five connected components. One of the 
components is self-conjugate and crosses the real axis with a double point at q ~ 
2.5328721401, from which there emerges a small horizontal segment running from 
2.5286467909 < q < 2.5370979311 (see Figure |a for detail). The other four connected 
components form two mutually conjugate pairs. One pair consists of arcs running from 
q « 0.0689480595 ± 0.9874383424 % to q « 1.6648104050 ± 2. 1404062947 i. The other 
pair exhibits T points at q « 2.039± 1.964 i, q w 2.332±1.638z and q « 2.478± 1.213 i, 
the latter of which is the end of a small bulb-like region enclosed by the limiting curve 
(see Figures |8|b,c). There are two small gaps at q w 1.67 ± 2.14 i. 

There are 14 endpoints: 

q w 0.0689480595 ±0.9874383424? (5.13a) 
q w 1.6648104050 ±2. 1404062947 i (5.13b) 
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g « 1.6870381566 ± 2. 1423501191 i (5.13c) 

q » 2.0370674106 ± 1.9742433636 i (5.13d) 

g » 2.3334923547 ± 1.6492963460 i (5.13e) 

g « 2.5286467909 (5.13f) 

g « 2.8373380200 ± 0.6533586125 i (5.13g) 

g w 2.5370979311 (5.13h) 

The bulb-like region is rather unusual. The point at which it starts, q ~ 2.478 ± 
1.213 i, really is a T point: computations show that each of the three curves of dom- 
inant equimodularity arriving at the T point has a smooth continuation into the 
region beyond the T point, where it becomes a curve of subdominant equimodular- 
ity (see Figure ^|c, where the dominant curves are shown with solid black lines and 
the subdominant curves with dashed red lines); moreover, the t value varies contin- 
uously along each of these three curves. One of the subdominant curves lies entirely 
within the enclosed bulb-like region and ends at a t = subdominant endpoint 
q ~ 2.5018915620 ± 1.15014645061 This means that there is a branch cut for these 
two subdominant eigenvalues, which become the dominant eigenvalues outside the 
bulb-like region. 

By Lemma |2]2|, the points where at least one amplitude vanishes are given by the 
condition det-D(g) = 0. This determinant can be written as 

detD(q) = g 13 (g- l) 23 (g- 2) 12 (g 2 -3g+ l) 8 (g- 3) 27 (g 3 -5g 2 + 6g- l)P(g) 2 (5.14) 

where P(q) is a polynomial of degree 218 with integer coefficients that we report in 
the file transfer 1 .m. There are two trivial zeros q — 0, 1 where all the amplitudes 
vanish simultaneously. There are five non-trivial zeros of det D(q) that correspond to 
the vanishing of the dominant amplitude: q = 2, q w 2.0617791396 ± 1.7315562279 i, 
q « 2.3406021969 ± 1.3825644365 % (shown with a x in Figure fl). 

The first non-trivial real zero converges rapidly to the Beraha number = 2 
(see Table |3|), while the second non-trivial zero appears to be converging (at a 
roughly 1/n 2 rate) to the real endpoint q ~ 2.5286467909. We expect that there 
will be further real zeros (whose number increases with n) that tend to the segment 
[2.5286467909 . . . , 2.5370979311 . . .] of the limiting curve and in the limit n -> oo 
become dense on that segment. The convergence to the complex isolated limiting 
points at q « 2.0617791396 ± 1.7315562279? and q « 2.3406021969 ± 1.3825644365 i 
is in both cases very fast and is compatible with an exponential rate. 

It is curious that (|5.14 ) vanishes also at the Beraha numbers B 5 = (3+v^5)/2, Bq = 



3 and B-j w 3.246979603717 — and hence also at their conjug ates B* = (3 - VE)/2, 
Bf ] w 1.5549581321 and Bf ] w 0.1980622642 — though all of these correspond to 
the vanishing of a subdominant amplitude. 



5.6 L x = 7f 

The transfer matrix is 32-dimensional; it can be found in the Mathematica file 
transfer 1 .m. 
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The limiting curve B (see Figure ||) appears to have seven connected components 
(but there may be more: we cannot be sure, as we were unable to compute the t = 
resultant). One of these components is a self-conjugate arc, and crosses the real axis 
at g ~ 2.6062482130; it has endpoints at g « 2.622974 ± 0.6095481 The other six 
components form three mutually conjugate pairs. One pair consists of arcs running 
from q « -0.002412 ± 0.933080 i to q « 1.425603 ± 2.248902 i. A very small gap (see 
Figure |TD|a for detail) separates these arcs from the second pair, which starts at the 
endpoint q « 1.433184 ± 2.248834 i and ends in a tiny bulb-like region with a T point 
at q « 2.415 ± 1.497 % (see Figure |l0|b). A small gap (see Figure |l0lc) separates these 



components from the third pair, which exhibits T points at q ~ 2.577 ± 1.133 i and 
endpoints at q ~ 2.616006 ± 0.616910 i; the latter is in turn separated by a very small 
gap from the self-conjugate arc (see Figure |10|d) . 

Please note that B enters for the first time into the half-plane Reg < 0; we 
conjecture that this occurs for all lattice widths L x > 7p. 

There are probably 14 endpoints: 



q 




-0.002412 ± 0.933080 i 


(5.15a) 


Q 




1.425603 ± 2.248902 i 


(5.15b) 


Q 




1.433184 ± 2.248834 % 


(5.15c) 


Q 




2.445207 ± 1.471332 % 


(5.15d) 


Q 




2.616006 ±0.616910i 


(5.15e) 


1 




2.622974 ± 0.609548 i 


(5.15f) 


Q 




2.886041 ± 0.602908 i 


(5.15g) 



We have determined these endpoints by the direct-search method; therefore, they 
are less accurate than endpoints computed by the resultant method, and the list is 
possibly incomplete. 

Due to limitations of CPU time, we were unable to obtain the explicit expression 
for det-D(g) as a polynomial in g; we were therefore unable to obtain all its roots. 
Instead, we evaluated det D(q) numerically at selected values of g. In particular, mo- 
tivated by the results of the preceding subsections, we computed det D(q) numerically 
at the first 50 Beraha numbers; and in those cases where det D(q) = 0, we numeri- 
cally diagonalized the transfer matrix in order to ascertain which amplitude (s) are the 
one(s) that vanish. We find two trivial zeros q = 0, 1 where all the amplitudes vanish 
simultaneously, and one non-trivial real zero q = 2 where the dominant amplitude 
vanishes (along with others). In addition, det D(q) vanishes at the Beraha numbers 
B 5 = (3 + y/E)/2, B 6 = 3, B 7 w 3.246979603717 and B 8 = 2 + ^2 (and hence also at 
their conjugates); these all correspond to the vanishing of a subdominant amplitude. 
Finally, inspection of Figure |9| suggests the existence of at least two pairs of complex- 
conjugate isolated limiting points, q ~ 2.48873±0.75416 i and q ~ 1.65436±2. 01881 i; 
and we confirm that the absolute value of the dominant amplitude is very small in 
both cases (2.9 x 10~ n and 2.1 x 10~ 6 , respectively), suggesting that this amplitude 
does indeed have a zero nearby. There might exist further isolated limiting points 
not found here. 
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The first non-trivial real zero converges quickly to the Beraha number B4 = 2 
(see Table §), and the next real zero converges (at a roughly 1/n rate) to the value of 
go ~ 2.6062482130 for this lattice, which is slightly smaller than the Beraha number 
B 5 . 



5.7 L x — 8p 

The transfer matrix is 70- dimensional; it can be found in the Mathematica file 
transfer 1 .m. 

The limiting curve B (see Figure |ll]) appears to consist of six connected compo- 
nents that form three mutually conjugate pairs (but there may be more). One pair is 
defined by arcs running from q w -0.054426±0.884363z to q w 1.211959±2.301760i. 
A very small gap (see Figure [jjja.) separates these arcs from the second pair of arcs, 
which run from q w 1.214531 ± 2.301385z to q w 2.326565 ± 1.753667 i. Another 
small gap (see Figure |l2"lb) separates these arcs from the third pair of components, 
which run from q w 2.330755 ± 1.737504 i to q w 2.660260 ± 0.001257 i and also have 
T points at q ~ 2.640±1.114i (see Figure |l2]c) . Note that the limiting curve does not 
cross the real axis: the closest approach is q ~ 2.660260 ± 0.001257 % (see Figure |12]d) . 

There are probably 14 endpoints: 



q « 


-0.054426 ± 0.884363 z 


(5.16a) 


g ps 


1.211959 ± 2.301760 1 


(5.16b) 


g w 


1.214531 ± 2.301385 z 


(5.16c) 


g ps 


2.326565 ± 1.753667 z 


(5.16d) 


g ps 


2.330755 ± 1.737504 i 


(5.16e) 


g « 


2.660260 ± 0.001257 i 


(5.16f) 


g ^ 


2.921658 ± 0.560969 i 


(5.16g) 



We have again determined these endpoints by the direct-search method. 

Once again we were unable to compute the determinant det D(q) as a polynomial 
in q, so we followed the same numerical method as in Section Exolto locate at least some 



of the isolated limiting points. There are of course two trivial isolated real zeros at q = 
0, 1 where all the amplitudes vanish simultaneously. In addition, det D(q) vanishes at 
the Beraha numbers B 4 = 2, B 5 = (3 + VE)/2, B 6 = 3, B 7 « 3.246979603717, B 8 = 
2 + V2 and B 9 « 3.5320888862 (and hence also at their conjugates). Unfortunately, 
we were unable to compute the corresponding amplitudes: even using 4000-digit 
numerical precision, the inverse of the change-of-basis matrix was obtained with no 
precision at all! Nevertheless, by inspection of Figure || we can guess that B4 and 
£> 5 are indeed isolated limiting points (see also the discussion below about the rate 
of convergence to those points), while B 6 , B 7 , B 8 , B 9 are not. Inspection of Figure |9] 
also suggests that there are two possible pairs of complex-conjugate isolated limiting 
points, namely q w 1.33321 ± 2.15164 i and q w 2.39242 ± 1.11180 i; but we are again 
unable to confirm this by a direct evaluation of the amplitudes at those points. 

The convergence of the first non-trivial real zero to B^ = 2 is very rapid (see 
Table H). The second non-trivial real zero converges fairly rapidly to the Beraha 
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number B$ = (3 + \/5)/2, which lies slightly below the point q ~ 2.66 where the 
limiting curve comes close to (but does not cross) the real axis. This rapid convergence 
suggests that B^ and B 5 are indeed isolated limiting points. This is the first square- 
lattice strip with free b.c. for which B$ appears as an isolated limiting point. 

6 Numerical Results for the Square-Lattice Chro- 
matic Polynomial: 

Periodic Transverse Boundary Conditions 

We have also computed the transfer matrix T(q) and the limiting curves B for 
square-lattice strips of widths 2 < L x < 8 with periodic boundary conditions in 
the transverse direction. We have checked our results for widths 2 < mp < 8 and 
lengths np = 2, 3 by comparing to the results of Biggs-Damerell-Sands ||109|| (resp. 



Shrock-Tsai f24] , f2~9[j ) for width np = 2 (resp. np = 3) and length mp, using the trivial 
identity 

Z(mp x np) = Z(np x mp) . (6-1) 

6.1 L x = 2p 

This is of course identical to L x = 2 F (Section phT| ). 



6.2 L x — 3p 

This case is also trivial, as the transfer matrix is one-dimensional. The result is 

Z(3 P x n F ) = q(q - l)(q - 2)(g 3 - 6q 2 + Uq - 13) n ~ 1 (6.2) 

The dominant-eigenvalue-crossing curve is of course the empty set B = 0. However, 
there are zeros for all n at q = 0, 1, 2 (trivially) and for all n > 2 at q ~ 1.7733011742± 
1.4677115087 i and q w 2.4533976515. 



6.3 L x = 4p 

The transfer matrix is two-dimensional. In the basis P = {1, 813 + 624} it can be 
written as 

rp/ a \ ( g 4 -8g 3 + 28g 2 -51g + 41 2(g 3 - 6g 2 + 14g - 12) \ 

T(4p) " { 2q-5 g 2 -4g + 5 J ' (6 ' 3) 



and the partition function is equal to 



g 2 -3g + 3\ T , m _ x ( 1 



Z(4 P xn F ) = q{q-l) ( ) • T(4p)™- 1 • ( J ) . (6.4) 
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The eigenvalues of the transfer matrix T(4p) and their corresponding amplitudes 
are given by (^13D / (^14D /(^20|) with 31 



P x (q) = q 4 - 8g 3 + 29g 2 - 55g + 46 (6.5a) 
p 2 ( q ) = g 8 -16g 7 + 118g 6 - 526g 5 + 1569g 4 - 3250g 3 + 4617g 2 

-4136g+1776 (6.5b) 

P 3 (q) = g(g-l)(g 6 -llg 5 + 54g 4 -152g 3 + 266g 2 - 277g + 128) (6.5c) 

P A {q) = 2g 2 (g-l) 2 (g-2)(g 2 -3g + l)(2g-5) 2 (6.5d) 

The limiting curve B (see Figure |T^) contains three pieces. Q It crosses the real 
axis at the double point q « 2.3026282864 (with t w 0.092), from which there emerges 
a horizontal segment running from g « 2.2533697671 to q ~ 2.3516882809. 

There are eight endpoints given by the zeros of the resultant at t — [= — P2(q)]: 

q « 0.7098031013 ± 2.0427103451 i (6.6a) 

g « 1.9923366166 ± 1.5941556425 i (6.6b) 

q « 2.2533697671 (6.6c) 

g « 2.3516882809 (6.6d) 

g « 2.9953312581 ± 1.4266372190 i (6.6e) 



The zeros of the amplitudes can be found by solving P±(q) = 0. There are trivial 
zeros q = 0,1, where both amplitudes vanish simultaneously; both these zeros lie 
in the region where the eigenvalue A + is dominant. There are five non-trivial zeros: 
q = 2, 5/2 (double) and (3± v / 5)/2. All of them lie in regions where there is a unique 
dominant eigenvalue; but only for q = 2 does the amplitude corresponding to the 
dominant eigenvalue vanish. The Beraha-Kahane-Weiss theorem thus implies that 
the isolated limiting points of zeros are precisely q = 0, 1,2. From Table ^ we see 
that the first non-trivial real zero does indeed converge rapidly to the Beraha number 

= 2. In addition, there are further real zeros (whose number increases with n) 
that tend to the segment [2.2533697671 . . . , 2.3516882809 . . .] of the limiting curve 
and in the limit n — ► 00 become dense on that segment. 

It is again curious that the amplitude corresponding to the subdominant eigen- 
value vanishes at the Beraha number = (3 + v / 5)/2 and its conjugate = 
(3-V5)/2. 

6.4 L x = 5p 

The transfer matrix is again two-dimensional. In the basis P 
£35 + #4i + <$52 } it can be written as 

rgp) = 

31 Our results are in agreement with those of eqns. (7.1)-(7.4)]. In particular, we have 
&sg(4), c ,i = -trT and 6 sg (4), c ,2 = detT. 

32 This curve is also depicted in jl9|, Figure 3(a)]. 



= {1, 5l3 +<$24 + 
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g 5 - 10g 4 + 45g 3 - 115g 2 + 169g - 116 5(g 4 - 9g 3 + 34g 2 - 63g + 47) 
g 2 - 6g + 10 g 3 - 9q 2 + 29q - 32 



(6.7) 



and the partition function is equal to 



Z(5 P x m F ) = q(q - l)(q - 2) ( + 2 ^ • T(5 P ) m - 1 • ( J ) . (6.8) 

The eigenvalues of the transfer matrix T(5p) and their corresponding amplitudes 
are given by (glg)/(glg)/(ggg) withQ 

P x ( g ) = g 5 - 10g 4 + 46g 3 - 124g 2 + 198g - 148 (6.9a) 
P 2 ( g ) = g 10 -20g 9 + 188g 8 - 1092g 7 + 4356g 6 - 12596g 5 + 27196g 4 

-44212g 3 + 52708g 2 - 41760g + 16456 (6.9b) 
P 3 (q) = g(g-l)(g-2)(g 7 -12g 6 + 66g 5 -214g 4 + 450g 3 

-646g 2 + 608g - 268) (6.9c) 

p 4 ( g ) = 5 g 2( ? _i)2( g _2) 2 (g 2 -3g+l)(g-3)(g 2 -6g + 10) 2 (6.9d) 

The limiting curve contains five pieces, and it crosses the real axis at go ~ 
2.6916837012.0 There are ten endpoints given by the zeros of the resultant at t = 

h -P a (g)]: 

g » 0.1650212134 ± 1.9190897717 i (6.10a) 

g « 2.0895893895 ± 1.9436539472 i (6.10b) 

g w 2.5034648023 ± 2.0851731765 i (6.10c) 

g w 2.5680063227 ± 0.4886738235 i (6.10d) 

g w 2.6739182721 ± 0.5983324603 i (6.10e) 

The zeros of the amplitudes can be found by solving P±(q) = 0. There are trivial 
zeros at q = 0, 1, 2, where both amplitudes vanish simultaneously; all of them lie in 
the region where the eigenvalue A_ is dominant. The non-trivial zeros of P4 are q = 3, 
3 ± i (each of them double) and (3 ± y/5)/2. Only for q = (3 + v^5)/2 = B 5 does the 
amplitude corresponding to the leading eigenvalue vanish. So the isolated zeros are 
expected to converge when the strip length goes to infinity to the first four Beraha 
numbers B 3 , . . . , B 5 = 0, 1, 2, (3 + VE)/2. 

From Table ^ we see that the third real zero is trivially equal to P4 = 2 (a 
cylindrical square lattice of odd width is not 2-colorable). The first non-trivial zero 
converges rapidly to the Beraha number B 5 , while the last real zero seems to be 
converging slowly (at a roughly 1/n rate) to the value g ~ 2.6916837012 where the 
limiting curve B crosses the real axis. 

The next Beraha number B§ = 3 is also a zero of P^q), but this corresponds to 
the vanishing of the subdominant amplitude. 

33 Our results are in agreement with those of j?2[ eqns. (3.8)-(3.10)]. 

34 This curve is also depicted in [B21 Figure 2]. 
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6.5 L x = 6p 

The transfer matrix is five- dimensional; it can be found in the Mathematica file 
transfer 1 .m. 



The limiting curve B (see Figure [T5p has three connected components. Two of 
them form a pair of mutually conjugate arcs, which extend for the first time into 
the half-plane Reg < 0. The third component is self-conjugate and crosses the real 
axis with a double point at q ~ 2.6110856839, from which there emerges a small 
horizontal segment running from 2.6089429411 < q < 2.6132283584 (see Figure [TB|a 



for detail). Note that this crossing lies slightly below the Beraha number B§ = 
(3 + y/E)/2m 2.6180339887. This component also has T points at q « 2.650 ± 1.240 i 
(see Figures [L6|b,c for detail). 

Inspection of Figure [15| might lead one to think that the self-conjugate component 
also exhibits a cusp at q « 2.568 ± 1.395 i: indeed, the curve B seems to have a dis- 
continuous derivative there. Moreover, there is no nearby dominant endpoint (which 
rules out the alternative hypothesis of a T point from which there emerges a very 
short curve terminating at a dominant endpoint); and anyway the value of t is also 
reasonably continuous around the alleged cusp (providing further evidence against 
the idea of a T point). However, by zooming on this region we can see that it is 
actually a single smooth curve (see Figure |I~6|b,c,d). What happens is that there is a 
subdominant endpoint very close to this curve (the subdominant curve is shown with 
dashed red lines); moreover, the subdominant eigenvalues are very close in modulus 
to the dominant ones. In particular, at the subdominant endpoint we have 

A dom = 10.4372144110 exp(-0.9719979634i) (6.11a) 
A sub = 10.4046337888 exp(-2.8500677201i) (6.11b) 

Near the subdominant endpoint, the two subdominant eigenvalues are very rapidly 
changing (as they have a square-root branch point), so that their crossing in modulus 
with the dominant eigenvalue occurs on a curve B of rapidly changing slope. 
There are ten endpoints: 

q « -0.1318891429 ± 1.7132242811 1 (6.12a) 

q w 1.9257517021 ± 2.2876287010 i (6.12b) 

q w 2.0571168133 ± 2.3885607275 i (6.12c) 

q « 2.6089429411 (6.12d) 

q « 2.6132283584 (6.12e) 

q « 3.1711921718 ± 0.8639071723 i (6.12f) 

By Lemma the points where at least one amplitude vanishes are given by the 
condition det D(q) = 0. This determinant can be written as 

det.D(g) = q 5 (q - l) 5 (g - 2) 4 (g 2 - 3q + l) 2 (q - 3) 3 (g 3 - 5q 2 + Qq - l)P(qf (6.13) 

where P(q) is a polynomial of degree 28 with integer coefficients that we report in 
the file transf erl .m. We find two trivial zeros q = 0,1 where all the amplitudes 
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vanish simultaneously. There are three points where a dominant amplitude vanishes: 
q = 2 and q « 2.4813444277 ± 1.7147613188 % (shown with a x in Figure [T|). The 
complex-conjugate pair of isolated limiting zeros lies extremely near, but not on, the 
limiting curve B (namely, about 0.006 to the left of B). The rate of convergence to the 
complex isolated limiting points is roughly 1/n instead of the expected exponential 
rate. This may be due to the fact that they are extremely close to the limiting curve 
B. 

From Table [| we see that the first non-trivial real zero converges rapidly to the Be- 
raha number i? 4 = 2, while the second non-trivial zero appears to be converging slowly 
(at a rate somewhere between 1/n and 1/n 2 ) to the real endpoint q « 2.6089429411. 
We expect that there will be further real zeros (whose number increases with n) that 
tend to the segment [2.6089429411 . . . , 2.6132283584 . . .] of the limiting curve and in 
the limit n — > oo become dense on that segment. 

It is curious that ( |6.13D vanishes also at the Beraha numbers B 5 = (3 + v / 5)/2, 
Bq = 3 and B 7 « 3.246979603717 and their conjugates, even though all of these 
correspond to the vanishing of a subdominant amplitude. 

6.6 L x = 7p 

The transfer matrix is six-dimensional; it can be found in the MATHEMATICA file 
transfer 1 .m. 

The limiting curve B (see Figure [17]) has seven connected components. One of 
them is a self-conjugate arc that crosses the real axis at go ~ 2.7883775115, which 
lies for the first time above the Beraha number B<j = (3 + y/E)/2 ~ 2.6180339887; this 
arc has endpoints at q « 2.7618995071 ± 0.4693560083 i. The other six components 
form three pairs of mutually conjugate components. The first pair are arcs running 
from q « -0.2962497164 ± 1.5256077564 i to q « 1.6542262925 ± 2.4866235231 i, and 
thus have support on Reg < 0. The second pair are also arcs, running from q w 
1.6947007027 ± 2.5327609879 i to q » 2.6589962013 ± 1.5245516751 i. The last pair of 
components has endpoints at q « 2.7275004011 ± 1.4172937300 i, q « 2.7873170476 ± 
0.4754613769 i and q w 2.8390155832 ± 1.3872842928 i, and T points at q w 2.737 ± 
1.405 i (see Figure |18|a,b); it is separated from the self-conjugate arc by a very small 
gap (see Figure |i~8|c). 

There are 16 endpoints: 



q 




-0.2962497164 ± 1.5256077564 i 


(6.14a) 


q 




1.6542262925 ± 2.4866235231 i 


(6.14b) 


q 




1.6947007027 ± 2.5327609879 i 


(6.14c) 


q 




2.6589962013 ± 1.5245516751 i 


(6.14d) 


q 




2.7275004011 ± 1.4172937300 i 


(6.14e) 


q 




2.7618995071 ± 0.4693560083 i 


(6.14f) 


q 




2.7873170476 ± 0.4754613769 i 


(6.14g) 


q 




2.8390155832 ± 1.3872842928 i 


(6.14h) 
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By Lemma |2]2|, the points where at least one amplitude vanishes are given by the 
condition det D(q) = 0. This determinant can be written as 



detD(g) = g 6 (g - l) 6 (g - 2) 6 (g 2 - 3g + l) 4 (g - 3) 3 (g 3 - 5g 2 + 6q - 1) x 

(q 2 - 4g + 2)P(q) 2 (6.15) 

where P(q) is a polynomial of degree 56 with integer coefficients. This polynomial can 
be found in the file transf erl .m. The trivial isolated limiting points are q — 0, 1, 2, 
where all the amplitudes vanish simultaneously. The non-trivial isolated limiting 
points are the Beraha number B$ = (3 + v^5)/2 and the pair of complex-conjugate 
roots q w 2.1027473746 ± 2.2083820861 i. This is the first square-lattice strip with 
cylindrical b.c. for which B 5 appears as an isolated limiting point. 

There is trivially a real zero at q = 2 (because the strip width is odd). The 
first non-trivial real zero (see Table ^) converges rapidly to the Beraha number B§ = 
(3 + \/5)/2, while the next real zero appears to converge slowly (at a roughly 1/n 
rate) to the value go ~ 2.7883775115 where the limiting curve B crosses the real axis. 
The rate of convergence to the complex isolated limiting points q ~ 2.1027473746 ± 
2.2083820861 i is very fast and is compatible with an exponential rate. 



It is curious that ( 6.15 ) vanishes also at the Beraha numbers B 6 = 3, B 7 m 
3.246979603717 and B 8 = 2 + ^2 and their conjugates, even though all of these 
correspond to the vanishing of a subdominant amplitude. 

6.7 L x = 8p 

The transfer matrix is 14- dimensional; it can be found in the Mathematica file 
transf erl .m. 

The limiting curve B (see Figure 19) has six connected components, which define 



three pairs of mutually conjugate components. The first pair is defined by arcs run- 
ning from q » -0.3908638747± 1.3698634697 i to q w 1.3863697070 ± 2.5801346584 i, 
which thus have support on Reg < 0. A small gap (see Figure pO|a) separates 
these components from the second pair, which consists of arcs running from q m 
1.3989312933±2. 5988401222 i to q w 2.5297861557±1. 8426263238 i. The last pair has 
endpoints at q w 2.5810431815 ± 1.8106192070 i, q w 2.7515311636 ± 0.0025313231 i, 
q w 2.7812812528 ± 1.0876657311 i and q w 3.2111321566 ± 0.6498638896 i, and T 
points at q w 2.801 ± 1.043 % and q w 2.783 ± 1.088 i. Note that the T points at 
q m 2.783 ± 1.088 i might look like cusps if we fail to use enough magnification (see 
Figure |20|b). But if we magnify the region sufficiently, we observe that it is indeed an 
ordinary T point (see Figure |20|c). Note also that B does not cross the real axis at 
any point; rather, there is a very tiny gap between q w 2.7515311636 ±0.0025313231 i 
(see Figure |20|d) . 

There are 16 endpoints: 

q w -0.3908638747 ± 1.3698634697 % (6.16a) 
q w 1.3863697070 ± 2.5801346584 i (6.16b) 
q w 1.3989312933 ± 2.5988401222 i (6.16c) 
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q « 2.5297861557 ± 1.8426263238 i (6.16d) 

q « 2.5810431815 ± 1.8106192070 i (6.16e) 

g « 2.7515311636 ± 0.0025313231 i (6.16f) 

q « 2.7812812528 ± 1.0876657311 i (6.16g) 

g « 3.2111321566 ± 0.6498638896 1 (6.16h) 

By Lemma |2]2|, the points where at least one amplitude vanishes are given by the 
condition det D(q) = 0. This determinant is given by 

det D(q) = q 14 (q-l) 20 (q-2) 13 (q 2 -3q + q 2 ) 11 (q-3) 18 (q 3 -5q 2 + Qq-l) 4 x 

(q 2 - 4g + 2)(g 3 - 6q 2 + 9q - l)P(q) 2 (6.17) 

where P(q) is a polynomial of degree 396 with integer coefficients (see file transf erl . m) 
The trivial isolated limiting points are q — 0,1: at these points all the ampli- 
tudes vanish simultaneously. The non-trivial isolated limiting points are the Be- 
raha numbers = 2 and B§ = (3 + \/5)/2 and the complex-conjugate pairs 
q « 1.6836371202 ± 2.4533856271 i and q w 2.6775096551 ± 1.2084144891 % (see Fig- 
ure [19]) . 

The first non-trivial real zero converges rapidly to the Beraha number B± = 2, 
and the next real zero converges rapidly to the Beraha number B 5 = (3 + y/E)/2 (see 
Table |j). The convergence to the complex isolated limiting points q ~ 1.6836371202± 
2.4533856271 i and g » 2.6775096551 ± 1.2084144891 i is also very fast and is com- 
patible with an exponential rate. 

It is curious that (|6.17| ) vanishes also at the Beraha numbers B 6 = 3, B 7 w 
3.246979603717, 5 8 = 2 + v 7 ^ and B 9 w 3.532088886238 and their conjugates, even 
though all of these correspond to the vanishing of a subdominant amplitude. 



7 Discussion and Open Questions 

7.1 Behavior of dominant-eigenvalue-crossing curves B 

In this paper we have computed the chromatic polynomials (= zero-temperature 
antiferromagnetic Potts- model partition functions) Pc{q) and their zeros for square- 
lattice strips of width 2 < L x < 8 and arbitrary length L y with free and cylindri- 
cal boundary conditions. In particular, we have extracted the limiting curves B of 
partition-function zeros when the length L y goes to infinity at fixed width L x . By 
studying the finite-width limiting curves and their behavior as we increase the width 
L x , we hope to shed light on the thermodynamic limit L x , L y — > oo. 

In Table [5] we summarize the main properties of the limiting curves B (and of the 
isolated limiting points) for all the lattices studied in the previous sections. Note the 
identity 

endpoints = (2 x components) + (2 x double points) + (T points) — 

(2 x enclosed regions) , (7.1) 
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which can be derived by simple topological/graph-theoretic arguments. 

Our first conclusion is that the limiting curves become in general more complicated 
as the strip width L x grows. In particular, the number of connected components, the 
number of endpoints, and the number of T points all tend to increase with the width 
L x . (Note that our counts for L x = 7p, 8f are only lower bounds on the true values.) 
Moreover, the size of the gaps between connected components, and the lengths of the 
protruding arcs associated to some of the T points, both seem to decrease with the 
strip width L x . The approach to the thermodynamic limit thus appears to be rather 
complicated. 

A second point deals with the existence or not of enclosed regions. Shrock p8 



Section III, point 3] conjectured that for families of graphs with a well-defined lattice 
structure, a sufficient condition for B to separate the g-plane into two or more regions 
is that the graphs contain at least one "global circuit" , defined as a route following 
a lattice direction which has the topology of S l and a length that goes to infinity as 
L y —>■ oo. (For strip graphs, this condition is equivalent to having periodic boundary 
conditions in the longitudinal direction.) Our results for L x = 6f, 7f show that this 
condition, whether or not it is in fact sufficient for the existence of enclosed regions, 
is in any case not necessary: enclosed regions can arise also with free longitudinal 
b.c.Q In both these cases, the enclosed regions are small bulb-like regions located at 
the end of one of the components of B. Enclosed regions do, however, seem to be 
atypical for square- lattice strips with free longitudinal b.c. 

A third point concerns the existence of chromatic zeros with Reg < 0, and more 
specifically of limiting curves B that intersect the half-plane Reg < 0. In 1980, Farrell 
|137|| conjectured, based on computations with small graphs, that all chromatic roots 



have Reg > 0. We now know that this conjecture is false |fL38[|; indeed, there are 



families of graphs whose chromatic roots, taken together, are dense in the whole 
complex g-plane ||40|| . Nevertheless, chromatic roots do seem to have a tendency to 
avoid the left half-plane, and it would be interesting to know why. Our computations 
show that for L x = 7p,8p and L x = 6p,7p,8p, the locus B does intersect the half- 
plane Reg < 0. Indeed, we conjecture that this happens for square-lattice strips of 
all widths L x > 7 (free b.c.) or L x > 6 (cylindrical b.c). 

Although the limiting curves get more complicated as the strip width L x grows, 
they do exhibit some regularities, as can be seen in Table [5] and Figures 21 and 22 



In Figure we superpose the limiting curves for all the square-lattice strips with 
free boundary conditions. In the leftmost part of the plot (Reg < 2), we see that the 
arcs behave monotonically: as the width L x increases, the corresponding arc moves 
outwards. In particular, minReg decreases monotonically with the width L x (see 
Table |5|). A similar behavior is observed in the right part of the plot for L x > 4: the 
limiting curves have similar shapes and they move monotonically to the right as L x 
grows. In particular, the point go where B crosses the real axis increases monotonically 
with the strip width, as does max Reg (see again Table [5]). In general, the shapes of 



the limiting curves look roughly similar to those obtained by Baxter [ 551 Figures 5 and 



35 For earlier examples showing that enclosed regions can arise also with free longitudinal b.c, see 
1, Figures 2(b), 3(b) and 4(a,b)]. 
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6] for the triangular lattice. We conjecture that the rightmost endpoints of B, which 
lie at q ~ 2.92 ± 0.56 i for L x = 8, will tend (slowly) to close up at the critical value 
q c = 3 as L x — > oo. We also conjecture, in analogy with the triangular lattice, that 
the values go will tend to a number strictly less than 3, probably somewhere around 
2.9. However, our strip widths are still too small to give unambiguous evidence for 
or against these conjectures. 

In Figure E2I we show the limiting curves for all the square-lattice strips with 



cylindrical boundary conditions. In the leftmost part of the plot (Reg < 2), the 
behavior of the arcs is again monotonic in the strip width; in particular, min Re q is 
again a decreasing function of L x (see Table [|). However, the qualitative behavior of 
the limiting curves on the right side of the plot (Reg > 2.3) is clearly not monotonic: 
there is a notorious difference between strips with even width and those with odd 
width. This difference is, in fact, to be expected: with periodic transverse boundary 
conditions, odd widths are in some sense "unnatural" as they introduce frustration in 
the antiferromagnetic Ising system (i.e. they make the chromatic number 3 rather than 
2, thereby forcing a chromatic zero at q = 2). It is curious that the difference between 
even and odd widths is significant only on the rightmost part of the limiting curve 
(namely, the part nearest q = 2 and q = 3). In any case, if we consider the even and 
odd subsequences separately, then g and max Re g are again monotonically increasing 
functions of L x (see Table [|). Moreover, the limiting curves for the square-lattice 
strips with even width and cylindrical boundary conditions have, apart from the case 
L x = 4, a qualitative shape in agreement with that for free boundary conditions. 
Finally, let us compare the limiting curves for free and cylindrical boundary con- 



ditions (Figures ^T] and [22] ). If we focus on the leftmost part of the plots, we see 
that both sets of curves tend to larger values of |Img| as L x grows, but the curves 
for cylindrical boundary conditions reach large values of |Img| much faster. Likewise, 
on the rightmost part of the plots, both sets of curves tend to larger values of Re q 
as L x grows (modulo the even-odd oscillation for cylindrical b.c), but the curves for 
cylindrical boundary conditions do so somewhat faster. This suggests that the ther- 
modynamic limit is achieved faster with cylindrical boundary conditions than with 
free boundary conditions. 

Let us conclude by mentioning briefly the work of Bakaev and Kabanovich ||139|| , 
who computed the large-g series for the infinite-volume limiting chromatic polyno- 
mial through order z m [where z — l/(g — 1)], using the finite-lattice method. Tony 
Guttmann (private communication) has kindly analyzed their series using differ- 
ential approximants. He finds that the nearest singularity to z = lies at z ~ 
(0.155 ± 0.005) ± (0.37 ± 0.03)i, corresponding to q « (1.96 ± 0.15) ± (2.30 ± 0.15)i. 
Curiously enough, this corresponds quite closely to a gap (pair of nearby endpoints) 
for width 6p (see Section |6.5| ), where the limiting chromatic polynomial is indeed sin- 
gular. But this may be a coincidence, as there is no corresponding endpoint for width 
8p, and the nearby endpoints for strips with free b.c. do not seem to be converging to 
this value. ^] The large-g series shows no hint of singularity at q c = 3, which is very 

36 For 5f there is a gap around z « 0.20 ±0.41 i, but as the width increases, the real part decreases 
beyond the predicted singularity (0.13 for 6f, 0.08 for 7f, 0.04 for 8f). 
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likely a weak essential singularity. 

7.2 Behavior of amplitudes and the Beraha conjecture 

For all the lattices we studied (up to width L = 8), we observed empirically that 
there is at least one vanishing amplitude Qjj(g) at each of the Beraha numbers up to 
B L+1 . It is reasonable to conjecture that this holds for all L: 

Conjecture 7.1 For a square-lattice strip of width L with free or cylindrical bound- 
ary conditions, at each Beraha number q = B 2 , . . . , -Bl+i there is at least one vanishing 
amplitude a^q). That is, det D(q) = for q = B 2 , . . . , £>l+i- 

For all the cases we studied except 7f and 8f [where we were unable to compute 
an explicit expression for det D(q)], we also verified that none of the roots of det D(q) 
correspond to Beraha numbers beyond Bl+x- We conjecture that this holds for all L: 

Conjecture 7.2 For a square-lattice strip of width L with either free or cylindrical 
boundary conditions, det D(q) ^ for all q = B^ with k > L + 1. [We assume, 
of course, that the chromatic polynomial is written in such a way that there are no 
identically vanishing amplitudes.] 

In some cases we found real roots of det D(q) that are very close to Beraha numbers: 
for instance, the strip with L = 5p (resp. L = 6p) has a zero of det-D(g) very close 
to B-j (resp. B 9 and -829)- Moreover, det.D(g) does in general have many real roots 
with q > B L+ i (and indeed with q > 4), as well as real roots with q < 0. 

We can strengthen the preceding conjecture to assert that det D(q) is strictly 
positive at the Beraha numbers beyond Bl+i'- 

Conjecture 7.3 For a square- lattice strip of width L with either free or cylindrical 
boundary conditions, det D(q) > for all q = B^ with k > L + 1. 

We have verified this conjecture for all square-lattice strips with L <7 F and L < 8p up 
to k = 50. Furthermore, the function f(k) = det D(Bk) seems to be a monotonically 
increasing function of k for k > L + 1. 

Conjecture [7^] asserts that at each of the Beraha numbers B 2 , . . . , B L+1 there 
is at least one vanishing amplitude ati(q), but it says nothing about whether the 
vanishing amplitude belongs to a dominant or a subdominant eigenvalue. Basing 
ourselves on a suggestion of Baxter [^, p. 5255], we conjectured that at each Beraha 
number q — B 2 , . . . , -Bl+i, the amplitude a* (5) corresponding to the eigenvalue \*(q) 
vanishes: here \*(q) is the eigenvalue that is dominant at small real q (e.g. at q — 1), 
analytically continued up the real axis (or, if there is a branch point on the real 
axis, then just above or below the real axis). [Note that \*(q) remains dominant 
until the path of analytic continuation crosses £>; after that, it becomes in general 
subdominant.] We tested this conjecture numerically as follows: Choose a path in the 
complex g-plane starting at q = 1 and proceeding to the right just above or below the 
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real axis (note that for 8f and 8p we must keep |Img| < 0.001 in order to avoid going 
around an endpoint).^] Subdivide this path into very small steps, and at each point 
compute the eigenvalues of the transfer matrix, following A*(g) "by continuity" (that 
is, starting with A*(l), choose at each q value the eigenvalue that is closest to the 
one chosen at the preceding q value). Then, when the Beraha number Bk is reached, 
diagonalize the transfer matrix, rotate the left and right vectors u(q) and v(q), and 
test whether the amplitude a*(q) vanishes. 

We were unable to carry out this computation for L = 7 F and 8p: in the for- 
mer case, because of CPU-time limitations, and in the latter case, because of loss 
of numerical precision in diagonalizing the matrices (even when we used 1000-digit 
arithmetic). In all but one of the other cases (namely, for L < 6p and L < 7p), we 
found that a*(q) does indeed vanish (sometimes along with other amplitudes) at all 
the Beraha numbers up to B L+ i. The strips with L = 8p and 8p are, however, a 
different the limiting curve B does not cross the real axis. If we follow a path 

from q — 1 to q — +00 along the real axis, we do not cross B; so A*(g) stays dominant 
everywhere on the real axis (and in particular at the Beraha numbers Bq, . . . ,Bg). 
Therefore, if a*(g) were to vanish at B e , . . . , B$, as our conjecture asserts, then those 
Beraha numbers would be isolated limiting points. But Figures |TI| and IT5 show clearly 



that they are not! So our conjecture must be false for the strips 8p and 8p. (Indeed, 
for 8p we confirmed explicitly that a* vanishes at B%, . . . , B5 but not at Bq, . . . , Bg; 
rather, it is a subdominant amplitude that vanishes at the latter points.) More gen- 
erally, we can expect our conjecture to be false whenever B does not cross the real 
axis.Q We are therefore obliged to modify our conjecture as follows: 

Conjecture 7.4 For a square-lattice strip of width L with free or cylindrical bound- 
ary conditions, let A*(g) be the eigenvalue that is dominant at small real q (e.g. at 
q = 1), analytically continued up the real axis (or, if there is a branch point on the real 
axis, then just above or below the real axis). Then, provided that the limiting curve 
B crosses the real axis, at each Beraha number q = B2, . . . , the corresponding 

amplitude a*(g) vanishes. [Other amplitudes may vanish as well] 

All our numerical evidence is consistent with Conjecture |7.4| . 

Remark. In the case L = 8p, suppose that we follow a path that does cross the 
limiting curve B slightly above the endpoint at q « 2.7515 + 0.0025 i. Then the 
eigenvalue A* will cease to be dominant to the right of B\ and we find that the 

37 Any two such paths give the same analytic continuation provided that there are no endpoints 
in the region between them. 

38 For the same reason, if in the cases L < 6f and L < 7p we analytically continue the dominant 
eigenvalue at q = 1 (namely, A*) to the region q > qo(L) by following a path that does not cross B 

— e.g., by going through one of the gaps in B — then the analytic continuation of A» will remain 
dominant everywhere. But then, the analytic continuation of the amplitude a» will not vanish at 
any of the Beraha numbers to the right of B, since we know that these are not isolated limiting 
points. Rather, the analytic continuation of one or more sufrdominant amplitudes at q = 1 will 
vanish at those Beraha numbers. We have numerically tested this behavior of the eigenvalues and 
amplitudes along paths that do not cross B in the same way we did for the paths that do cross B. 
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corresponding amplitude a* does vanish at all the Beraha numbers Bq, . . . , Bg\ So a 
variant of Conjecture [T4] applies in this well. But we are unable to see what 

general principle might be at work. 

Conjecture [T4| "explains" why the first few Beraha numbers — but only the first 
few — arise as limiting points of chromatic roots, at least in those cases where the 
limiting curve B crosses the real axis at some point qo(L). Indeed, those Beraha 
numbers that satisfy both q < -Bl+i and q < qo(L) correspond to the vanishing of 
the dominant amplitude, hence are isolated limiting points of chromatic roots. By 
contrast, the remaining Beraha numbers correspond either to the vanishing of a sub- 
dominant amplitude (in case qo(L) < q < B L+1 ) or do not correspond to the vanishing 
of any amplitude (in case q > B L+ i, assuming the validity of Conjecture |7~2]) -PI As 
the strip width L grows, the limiting curve B moves to the right and "uncovers" more 
Beraha numbers; qo(L) presumably tends to a limiting value go(°°)- F° r the triangu- 
lar lattice, Baxter's analytic solution predicts that qo(oo) ~ 3.81967, which lies 
between 5 14 and B i5 and in particular lies strictly below the critical point q c = 4. For 
the square lattice, an analytic solution is lacking, but our results in Table [5] suggest 
(assuming monotonicity in L) that qo(oo) > 2.788, and analogy with the triangular 
lattice suggests that go(oo) < q c = 3. It follows that go(oo) lies between B 5 and Bq, 
so that the first four Beraha numbers B 2 , . . . , B 5 — but only these — will be isolated 
limiting points. 



7.3 Upper zero-free interval for bipartite planar graphs 

Let G be a loopless planar graph. Then it is not hard to prove that Pg(q) > for 
all integers q > 5;F] moreover, one of the most famous theorems of graph theory - 



the Four-Color Theorem |[L41| , |T42| , [H3] , |144| |145|| - asserts that Pc(q) > holds in 
fact for all integers q > 4. 

It is natural to ask whether these results can be extended from integer q to real q. 
The answer is yes, at least in part: Birkhoff and Lewis ||101|| proved in 1946 that if G 



is a loopless planar graph, then Pc{q) > for all real numbers q > 5.0 Furthermore, 
they conjectured that Pc(q) > also for 4 < q < 5; and while no one has yet found 
a proof, no one has found a counterexample either, so it seems plausible (in the light 
of the Four-Color Theorem) that the conjecture is true. 

39 In case qo(L) is itself a Beraha number < Bl + i (as happens e.g. for L = 3f), it corresponds 
to case (b) of the Beraha-Kahane- Weiss theorem, hence to a non- isolated limiting point. More 
generally, if B intersects the real axis in an interval [qo^—^L), qo^(L)], then all the points in this 
interval correspond to non-isolated limiting points, even if a dominant amplitude should happen to 
vanish. 

40 This is the Five-Color Theorem, which goes back to Heawood in 1890. For a proof, see e.g. 
[140, Theorem V.8, pp. 154-155]; or for an elegant alternate proof of an even stronger result, see 



140, Theorem V.12, pp. 161-163]. 



41 See also Woodall ]102| , Theorem 1] and Thomassen 1 103 , Theorem 3.1 ff.] for alternate proofs 
of a more general result . 
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Now some planar graphs can be colored with three or even two colors; their chro- 
matic polynomials Pg(q) are strictly positive for integers q > 3 or q > 2, respectively. 
Can these bounds can be extended to real ql That is, if G is a /c-colorable planar 



graph, do we have Pg(q) > f° r ah rea l Q > k? Woodall ||102| , p. 142] conjectured that 
the answer is yes. For k — 4, this is the conjecture of Birkhoff and Lewis mentioned 
above. For k = 3, however, Thomassen [ 103|, pp. 505-506] has shown that Woodall's 



conjecture is false: there exist 3-colorable planar graphs with real chromatic roots 
greater than 3.0 We can now show that Woodall's conjecture is false also for k = 2: 
there exist 2-colorable (i.e. bipartite) planar graphs with real chromatic roots greater 
than 2. For example, the 4p x 6f square lattice has chromatic roots at q ~ 2.009978 
and q ~ 2.168344; and the same pattern persists for larger lattices, with both free 
and periodic transverse boundary conditions (see Tables |3| and f|). Indeed, for the 
cases 8p X n-p and 8p x %, we see numerically (Tables || and ^) that there are real 
chromatic roots tending to B 5 = (3 + y/E)/2 ~ 2.618034 from below as n — » oo.0 
This leads us to modify Woodall's conjecture as follows: 

Conjecture 7.5 Let G be a bipartite planar graph. Then -Pg(<?) > for real q > 
£ 5 = (3 + V5)/2. 

Let us make two remarks: 

1) Planarity here is crucial, as non-planar bipartite graphs can have arbitrarily 
large real chromatic roots. Indeed, the complete bipartite graphs K ni)Tl2 , in the limit 
n 2 — > oo with ni fixed, have real chromatic roots arbitrarily close to all the integers 
from 2 through |_ n i/2j ||104j , Theorem 8]. 

2) Initially we conjectured that not only Pg(q) but also all its derivatives are 
positive for q > B>§. But this turns out to be false: for example, the 8p x 16f lattice 
has P G (B 5 ) < 0; and the 8 P x 24 F lattice has P' G {q) < for 2.638337 < q < 2.687058. 

42 Start with a graph K and a real number qo for which Pk (go) < 0. Then Thomassen flO 



Theorem 3.9] constructs a 2-degenerate (and hence 3-colorable 14C, Theorem V.l, p. 148]) graph 
K(m) such that PK(m)(lo) < 0; moreover, K(m) can be chosen to be planar if K is. Since there 
exist planar graphs K with real chromatic roots q\ greater than 3, and since the 3-colorability of 
K(m) [or alternatively the Four-Color Theorem] implies that -Pftr( m )(4) > 0, we can take qo = qi — e 
and conclude that K(m) has a chromatic root in the interval qo < q < 4. Thus, the upper zero-free 
interval for 3-colorable planar graphs is the same as that for all planar graphs. 



In presenting this result, Thomassen 1 103 , p. 506] further asserted that there exist planar graphs 
K with real chromatic roots arbitrarily close to 4 j but this assertion apparently arises from a 
misunderstanding of the Beraha-Kahane p8| theorem that 4p x n-p triangular lattices have complex 
chromatic roots arbitarily close to 4. In fact we do not know of any planar graphs with real chromatic 
roots arbitarily close to 4. In our study of triangular-lattice strips ]87j we have thus far found 
chromatic roots up to ps 3.51; and Baxter's [|55| result qo{oo) ps 3.81967 suggests that sufficiently wide 
and long pieces of the triangular lattice will have real chromatic roots up to at least w 3.801938. 
But we do not know of any planar graphs with chromatic roots in the interval [i?i4, oo). 

We thank Carsten Thomassen and Douglas Woodall for correspondence concerning these ques- 
tions. 

43 The 5p x np and 7p x tif lattices have chromatic roots larger than B5, but graphs of odd strip 
width with periodic b.c. are not bipartite. 
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7.4 Prospects for future work 



One very interesting extension of this work is the computation of the chromatic 
polynomials for strips with periodic boundary conditions in the longitudinal direction 
Such computations have been performed for strip widths m = 2, 3 by ad hoc 



methods |[L09| , |34| , |35| , |23| , [24] , |26| , |27| , [291 , but a systematic transfer-matrix formalism 
has heretofore been lacking. In fact, the needed formalism can be obtained by a slight 
extension of the methods explained in this paper. Suppose we want to obtain the 
chromatic polynomial for a square-lattice strip mp/p x n-p (i.e., either free or periodic 
transverse b.c, and periodic longitudinal b.c). The idea is simple: instead of just 
keeping track of the connectivities among the m sites on the current top row, we keep 
track of the connectivities among the 2m sites on the current top row and the bottom 
row. Let us call these sites 1,2, ... , m and 1', 2', . . . , m', respectively. Initially the 
top and bottom rows are identical. We then enlarge the lattice one site at a time, 
exactly as in Section |3.2| ; the join and detach operations act on the sites of the top 
row, with those of the bottom row simply "going along for the ride". At the end, 
when we have obtained a lattice with n + 1 rows, we identify the top and bottom 
rows. The partition function for periodic longitudinal boundary conditions can thus 
be written as 

Z Gr (q,{v e }) = u T (VH)"v cquiv , (7.2) 

where "equiv" denotes the partition { {1, 1'}, {2, 2'}, . . . {m, m'} } and u T is defined 
by 

U = U S n i J 2 2' • " - -imm' (7.3) 

where u T is defined in (|3.41|) . 

In future work in collaboration with Jesper-Lykke Jacobsen, we will extend the 



results of the present paper to wider strip widths |p6|, to the triangular lattice [37 



to periodic boundary conditions in the longitudinal direction |p8[ , and to nonzero 
temperature [|65|. 
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n 


B n (exact) 


fi n (num) 


Pn{q) 


Other B y n ' 


2 








q 




3 


1 


1 


q-1 




4 


2 


2 


q-2 




5 


(3 + VE)/2 


2.6180339887 


q 2 - 3g + 1 


(3 - VE)/2 


6 


3 


3 


9-3 




7 




3.2469796037 


g 3 - 5g 2 + 6g - 1 


0.1980622642 
1.5549581321 


8 


2 + V2 


3.4142135624 


g 2 - 4g + 2 


2-V2 


Q 




O.OoZUooooOZ 


g — og -+- yg — i 


U. 1ZUD14 1 0o4 

2.3472963553 


10 


(5 + x/5)/2 


3.6180339887 


q 2 — 5g + 5 


(5 - V5)/2 


11 




3.6825070657 


g 5 - 9g 4 + 28g 3 - 35g 2 + 15g - 1 


0.0810140528 
0.6902785321 
1.7153703235 

/.ooUooUUZDU 


12 


2 + 


3.7320508076 


g 2 - 4g + 1 


2 - 


13 




3.7709120513 


g 6 - llg 5 + 45g 4 - 84g 3 + 70g 2 
-21g+ 1 


0.0581163651 
0.5029785037 
1.2907902259 
2.2410733605 
3.1361294935 


14 




3.8019377358 


g 3 - 7g 2 + 14g - 7 


0.7530203963 
2.4450418679 


15 


(9 + ^+^30-6^/4 


3.8270909153 


g 4 - 9g 3 + 26g 2 - 24g + 1 


0.0437047985 
1.7909430735 
3.3382612127 


16 


2 + a/2 + y/2 


3.8477590650 


g 4 - 8g 3 + 20g 2 - 16g + 2 


0.1522409350 
1.2346331353 
2.7653668647 



Table 1: Beraha numbers B n = 4cos 2 (7r/n) and their minimal polynomials p n {<l)- 
For each n we give the exact expression of the Beraha number B n whenever it can 
be expressed in terms of square roots alone; its numerical value to 10 decimal places; 
the unique irreducible monic polynomial p n (q) with integer coefficients having B n 
as a root; and the other zeros of p n (q), which are the generalized Beraha numbers 
B^ = 4cos 2 (A;7r/n) with k relatively prime to n. 
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m 


B(m) 


C 


TriFree(m) 


SqFree(m) 


dfn 


TriCyl(m) 


SqCyl(m) 


1 


1 


1 


1 


1 


l 


1 


1 


2 


2 


2 


1 


1 


l 


1 


1 


3 


5 


5 


2 


2 


l 


1 


1 


4 


15 


14 


4 


3 


3 


2 


2 


5 


52 


42 


9 


7 


6 


2 


2 


6 


203 


132 


21 


13 


15 


5 


5 


7 


877 


429 


51 


32 


36 


6 


6 


8 


4140 


1430 


127 


70 


91 


15 


14 


9 


21147 


4862 


323 


179 


232 


28 


22 


10 


115975 


16796 


835 


435 


603 


67 


51 


11 


678570 


58786 


2188 


1142 


1585 


145 


95 


12 


4213597 


208012 


5798 


2947 


4213 


368 


232 


13 


27644437 


742900 


15511 


7889 


11298 


870 


499 


14 


190899322 


2674440 


41835 


21051 


30537 


2211 


1241 



Table 2: Dimension of the transfer matrix. For each strip width m we give the num- 
ber B{m) of all partitions, the number C m of non-crossing partitions, the number 
TriFree(m) = M m _ x of non-crossing non-nearest-neighbor partitions with free bound- 
ary conditions, and the number SqFree(m) of equivalence classes of non-crossing non- 
nearest-neighbor partitions modulo reflection in the center of the strip. We also give 
the number d m (= R m for m > 2) of non-crossing non-nearest-neighbor partitions 
with periodic boundary conditions, the number TriCyl(m) of equivalence classes of 
such partitions modulo translations, and the number SqCyl(m) of equivalence classes 
of such partitions modulo translations and reflections. 
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Lattice 


3rd Zero 


4th Zero 


5th Zero 


6th Zero 


3 F X 3 F 


1.646039212420 








3 F x 6 F 










3 F x 9 F 


1.862295803794 








3 F x 12 F 










3 F x 15 F 


1.910244567418 








3 F x 18 F 










3 F x 21 F 


1.932253338339 








3 F x 24 F 










3 F x 27 F 


1.945103511556 








3 F x 30 F 










4 F X 4 F 










4 F x 8 F 










4 F x 12 F 


2.000607664038 


2.183434328589 






4 F x 16 F 


2.000017521546 


2.226186181588 






4 F x 20 F 


2.000000515361 


2.248253640526 






4 F X 24 F 


2.000000015170 


2.261494080470 






4 F x 28 F 


2.000000000447 


2.270172437566 






4 F x 32 F 


2.000000000013 


2.276213662199 






4 F x 36 F 


2.000000000000 


2.280609243979 






4 F x 40 F 


2.000000000000 


2.283918256290 






4 F x 100 F 


2.000000000000 


2.236070288638 


2.284202920228 


2.297805980307 


5 F x 5 F 


1.955073615801 








5 F x 10f 


2.000022457863 


2.243311545349 






5 F x 15 F 


1.999999994509 








5 F x 20 F 


2.000000000001 


2.335823711578 






5 F x 25 F 


2.000000000000 








5 F x 30 F 


2.000000000000 


2.365828458342 






5 F x 35 F 


2.000000000000 








5 F x 40 F 


2.000000000000 


2.381070502769 






5 F x 45 F 


2.000000000000 








5 F x 50 F 


2.000000000000 


2.390328275726 






6 F x 6 F 


2.001381451484 


2.196830038914 






6 F x 12 F 


2.000000000760 


2.390498998123 






6 F x 18 F 


2.000000000000 


2.448434501424 






6 F x 24 F 


2.000000000000 


2.475714120608 






6 F x 30 F 


2.000000000000 


2.491245543049 






6 F x 36 F 


2.000000000000 


2.501126630104 






6 F x 42 F 


2.000000000000 


2.507892809644 






6 F x 48 F 


2.000000000000 


2.512775536891 






6 F x 54 F 


2.000000000000 


2.516440469487 






6 F x 60 F 


2.000000000000 


2.519276871603 






6 F x 240 F 


2.000000000000 


2.534921463459 






7 F x 7 F 


1.999994176430 








7 F x 14 F 


2.000000000000 


2.451966225086 






7 F x 21 F 


2.000000000000 








7 F x 28 F 


2.000000000000 


2.523291736983 






7 F x 35 F 


2.000000000000 








7 F x 42 F 


2.000000000000 


2.548814353555 






7 F x 49 F 


2.000000000000 








7 F x 56 F 


2.000000000000 


2.562226841180 






7 F x 63 F 


2.000000000000 








7 F x 70 F 


2.000000000000 


2.570504933475 






8 F x 8 F 


2.000000005426 


2.391719919086 






8 F x 16 F 


2.000000000000 


2.531414423190 






8 F x 24 F 


2.000000000000 


2.576747844784 






8 F x 32 F 


2.000000000000 


2.597790566370 






8 F x 40 F 


2.000000000000 


2.609001821476 






o F x <±op 


z.uuuuuuuuuuuu 








8 F x 56 F 


2.000000000000 


2.617322619879 






8 F x 64 F 


2.000000000000 


2.617921204353 






8 F x 72 F 


2.000000000000 


2.618018253506 






8 F x 80 F 


2.000000000000 


2.618031848556 






Beraha 


2 


2.618033988750 







Table 3: Real zeros of the chromatic polynomials of finite square-lattice strips with 
free boundary conditions in both directions, to 12 decimal places. A blank means 
that the zero in question is absent. The first two real zeros q — 0, 1 are exact on all 
lattices. "Beraha" indicates the Beraha numbers B 4 = 2 and B 5 = (3 + \/5)/2. 
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Lattice 


3rd Zero 


4th Zero 


5th Zero 


6th Zero 


7th Zero 


8th Zero 


A \s A 
4p X 4p 














A \y Q 

4p X op 


z.UUUyo ( 04DDOO 


o oqq c^qoq t;i /I n/1 
z.Zoo0ozo0l4U4 










4p X 12p 


Z.UUUUllzyoool 


z.Zoolol24UlDy 










4p x 16f 


2. U0U0U0 139380 


o Qn'7t;ooooc;Q/iQ 
z.oU ( 0ZoZZ0o4o 










4p x 20p 


o nnnnnnnm vol 
2.U0U0U0U01 rzl 


o 91 rid Q^nonon 










4p x 24p 


2.UU0UUU0U0U21 


o o o t 1 o i o i r\ cr ~i ri 

z.oz74olol9olU 










4p x 28f 


2.UUUUUUUUUUUU 


o QQOt;QQnc;o/i'7i 
2.oo2oooUoo4/ 1 










4p x 32p 


2.U0U0U0U0U0U0 


o o o i o n o o i ri o o 










4p x 36f 


2.U0U0U0U0U0U0 


2.338792911735 










4p x 40f 


2.UUUUUUUUUUUU 


2 . <540o0785ood4 










^p X lUUp 


9 nnnnnnnnnnnn 
z.uuuuuuuuuuuu 




9 9708'}fi£89'3Qfi 
Z.Z lUOODDoZoyD 




9 Q(^1 1 QQQ97 


9 ^AQA9fi1 ^(^Q7S 


5p x 5f 


2 












5p x 10f 


2 


2.o7yDyz7y874d 










5p x 15f 


2 












5p x 20p 


2 


z.bloUoo / 4Zz4b 










5p x 25f 


2 


2.ol8482yyo587 


2. 6430458 14o23 








5p x 30p 


2 


z.di /yy4yyzzo4 










5p x 35f 


2 


2.618037byb771 


2.658yU8973824 








5p x 40f 


2 


2.618033639521 










5p x 45f 


2 


2.ol8Uo402lb7b 


2.6ob710728o80 








5p x 50f 


2 


2.0loUooyo004o 










Dp X OF 


Z.UUUUU44S4D ID 


z.4U / 4yooo /Uoz 










6p x 12p 


o nnnnnnnnnnnn 
2.UUUUUUUUUUUU 


o t;i<;c;if!inco/i7 
z.0lD0lbiyDZ4 / 










6p x 18f 


o nnnnnnnnnnnn 
2.UUUUUUUUUUUU 


o k c 1 /inccoconn/^ 

2.ooi4yooozyuo 










6p x 24p 


2.0U0U0U0U0U0U 


2.5o8o457104oo 










6p x 30f 


o nnnnnnnnnnnn 
2.UUUUUUUUUUUU 


o C7Q7/I Vi^nnnT'T' 










6p x 36f 


o nnnnnnnnnnnn 
2.UUUUUUUUUUUU 


z.OoooOoUoZDlo 










6p x 42p 


2.U0U0U0U0U0U0 


2.5yUUU8147yb5 










6p x 48p 


2.U0U0U0U0U0U0 


2.5yo4o5oooiyz 










6p x 54p 


o r\dr\dr\dr\dr\dr\d 

2.U0U0U0U0U0U0 


2.5yb0602ob523 










6p x 60f 


o r\r\r\dr\r\r\dr\dr\d 

2.0U0U0U0U0U0U 


o cnoi om ^1 cro7 

2.5y81zyiDloo7 










Dp X z4Uf 


o nnnnnnnnnnnn 

Z.UUUUUUUUUUUU 


9 fii n i 7cni^9"i son 










7p x 7f 


2 












7p x 14p 


2 


2.ol7y37723253 










7p x 21p 


2 


o £?1 ono 7707 

2.ol8Uo4U17737 


2.721810707015 








7p x 28f 


2 


O ftl 00990007/1 1 

2.oloUooyoo(41 










7p x 35f 


2 


o ono90oo7cri 

2.618033988750 


2.748ooz7dz81z 








7p x 42p 


2 


o ^1 onooooo^fn 

2.ol5Uooyo875U 










7p x 49f 


2 


o <;i oo99noo7c;o 


2. /bU2oUUoD0lo 








7p x 56f 


2 


o c \ onoonoo7rn 

2.ol8Uooyo875U 










7p x 63f 


2 


O 00990007KO 


o /inncrnonQcc 
2. /bb4yyoUoUoo 








i p A 1 uf 


2 


2 nl 80 ^9887 'Sf) 

iiUlOUUO JO O 1 tJL/ 










8 P x 8 F 


2.000000000001 


2.551072878420 










8 P x 16p 


2.000000000000 


2.616714700486 










8 P x 24 P 


2.000000000000 


2.618032009068 










8 P x 32 P 


2.000000000000 


2.618033986108 










8 P x 40 P 


2.000000000000 


2.618033988746 










8 P x 48 P 


2.000000000000 


2.618033988750 










8 P x 56 P 


2.000000000000 


2.618033988750 










8 P x 64 P 


2.000000000000 


2.618033988750 










8 P x 72 P 


2.000000000000 


2.618033988750 










8 P x 80 P 


2.000000000000 


2.618033988750 










Beraha 


2 


2.618033988750 











Table 4: Real zeros of the chromatic polynomials of finite square-lattice strips with 
periodic boundary conditions in the transverse direction and free boundary conditions 
in the longitudinal direction, to 12 decimal places. A blank means that the zero in 
question is absent. The first two real zeros q = 0, 1 are exact on all lattices; the third 
real zero q = 2 is exact on all lattices of odd width. "Beraha" indicates the Beraha 
numbers 5 4 = 2 and B 5 = (3 + v / 5)/2. 
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Eigenvalue-Crossing Curves B 


Isolated Points 


Lattice 


#C 


#E 


#T 


#D 


#ER 


min Re q 


qo 


max Re q 


#RI 


#CI 


2 F 























2 





3f 


3 


6 











0.586570 


2 


2.5 


2 





4f 


3 


10 


2 


1 





0.325474 


[2.228359,2.301416] 


2.667426 


3 





5f 


5 


14 


4 








0.170897 


2.428438 


2.769205 


3 





6f 


5 


14 


6 


1 


2 


0.063142 


[2.528647, 2.537098] 


2.837338 


3 


2 


7f 


7t 


14t 


4t 


ot 


2t 


-0.044443 


2.606248 


2.886041 


3t 


2t 


8f 


6t 


14+ 


2+ 


0+ 


ot 


-0.130642 
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Table 5: Summary of qualitative results for the eigenvalue-crossing curves B and 
for the isolated limiting points of zeros. For each square-lattice strip considered in 
this paper, we give the number of connected components of B (# C), the number 
of endpoints E), the number of T points T), the number of double points 
D), and the number of enclosed regions (# ER); we give the minimum value of Reg 
on B, the value(s) qo where B intersects the real axis (* denotes an almost-crossing), 
and the maximum value of Reg on B. We also report the number of real isolated 
limiting points of zeros (# RI) [which are always successive Beraha numbers B 2 , B 3 , 
. . .] and the number of complex-conjugate pairs of isolated limiting points (# CI). 
The symbol + indicates uncertain results. 
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Figure 1: Schematic representation of a generic T point where three eigenvalues 
Ai, A2, A3 are simultaneously dominant. The solid black lines represent the loci of 
dominant eigenvalue crossings, while the dashed red lines represent the loci of sub- 
dominant eigenvalue crossings. Each line has been labeled by the inequalities and 
equalities it satisfies, and each region has been labeled by the inequalities it satisfies. 
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Figure 2: Zeros of the partition function of the g-state Potts antiferro magnet on the 
square lattices 3^ x 15^ (squares), 3^ x 30f (circles) and 3f x oo^ (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Zeros sq lattice L = 4 
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Figure 3: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 4p x 20^ (squares), 4p x 40p (circles) and 4p x oo^ (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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(a) (b) 

Figure 4: Detail of the limiting curves B for the g-state Potts antiferromagnet on 
a square lattice 4 F x oop. (a) Region near the double point q ~ 2.2649418565. The 
value of t is continuous around the double point, with t « 0.0621. (b) Region near 
the T point q « 2.327 + 0.9113 i. At this T point we have t w (1.818, 12.962,0.655) 
and hence 9 w (2.136, 2.988, 1.160), so that £ 6 = 2vr. 
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Figure 5: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 5^ x 25^ (squares), 5^ x 50f (circles) and 5f x oo^ (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Figure 6: Detail of the limiting curves B for the g-state Potts antiferromagnet on 
a square lattice 5^ x oo^. Region near the T points of the limiting curve. On the 
upper T point q ~ 2.291 + 1.561 i, we have t ~ (0.999,0.179, 1.434), corresponding 
to 6 « (1.569,0.354,1.924). On the lower T point q w 2.423 + 0.1067?, we have 
t w (1.823, 0.434, 0.774) with 9 w (2.138, 0.820, 1.318). 
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Zeros sq lattice L = 6 
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Figure 7: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 6f x 30^ (squares), 6^ x 60f (circles) and 6f x oo^ (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Figure 8: Detail of the limiting curves B for the g-state Potts antiferromagnet on 
a square lattice 6^ x oop. (a) Region around the double point at q w 2.53287. The 
value of t is continuous around this double point, with t ~ 0.00985. (b) Region 
containing the three T points. On the upper T point q ~ 2.039 + 1.964 i, we have 
t w (0.871,0.0521,0.970) with 6 w (1.434,0.104,1.540); on the middle T point q w 
2.332 + 1.638 i, we have t w (18.021,0.0843,7.119) and 6 w (3.031,0.168,2.863); 
on the lower T point q w 2.478 + 1.213 i, we have t w (0.272,0.618,1.069) and 
9 w (0.532, 1.107, 1.638). (c) Detail of the bulb-like region around q w 2.478 + 1.2131 
Dominant crossing curves are depicted in solid black lines, while subdominant crossing 
curves are shown with dashed red lines. 
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Zeros sq lattice L = 7 
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Figure 9: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 7^x35^ (squares), If x 70f (circles) and 7f x oo^ (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the direct-search method. 
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Figure 10: Detail of the limiting curves B for the g-state Potts antiferromagnet on a 
square lattice 7f x oof- (a) Region around the gap between q w 1.425603 + 2.248902 i 
and q « 1.433184 + 2.248834 i. (b) Bulb-like region around the T point q m 2.415 + 
1.4971 At this point we have t « (0.023, 1.248, 1.310) and 6 « (0.047, 1.791, 1.838). 
(c) Region around the gap between the bulb-like region at q ~ 2.415 + 1.497 % and the 
endpoint at q w 2.445207+ 1.471332 z. There is also a T point at q w 2.577+ 1.133 z, 
where i w (2.108,2.737,1.016) and w (2.257,2.441,1.586). (d) Region around the 
gap between q w 2.616006 + 0.616910 1 and g w 2.622974 + 0.609548 i. 
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Zeros sq lattice L = 8 
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Figure 11: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 8^ x 40p (squares), 8^ x 80f (circles) and 8f x oo^ (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the direct-search method. 
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Figure 12: Detail of the limiting curves B for the g-state Potts antiferro magnet on 
a square lattice 8p x oof- (a) Region around the small gap at q ~ 1.21 + 2.30 i. 
(b) Region around the gap at q ~ 2.32 + 1.75 i. (c) Region around the T point 
at g k 2.640 + 1.114*. At this point we have t w (0.993,1.013,0.0113) and w 
(1.563, 1.584,0.023). (d) Region around the tiny gap at q ta 2.660260. 
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Zeros sq lattice L = 4 
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Figure 13: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 4p x 20p (squares), 4p x 40p (circles) and 4p x oop (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Zeros sq lattice L = 5 
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Figure 14: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 5p x 25p (squares), 5p x 50p (circles) and 5p x oop (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Zeros sq lattice L = 6 
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Figure 15: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 6p x 30^ (squares), 6p x 60f (circles) and 6p x oo? (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Figure 16: Detail of the limiting curves B for the g-state Potts antiferromagnet on a 
square lattice 6p x oop. (a) Region around the double point q ~ 2.6110857. At this 
double point we have t m 0.0053. (b) Region around the T point at q ~ 2.650 + 1.240 i. 
At this T point we have t « (0.392,1.224,0.562) and 6 w (0.748,1.771,1.024). (c) 
The same as in (b), but we show the dominant (solid black line) and subdominant 
(dashed red line) crossing curves, (d) Blow-up of region around the quasi-cusp at 
q « 2.568 + 1.3981 



84 



Zeros sq lattice L = 7 
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Figure 17: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 7p x 35^ (squares), 7p x 70 F (circles) and 7p x oo? (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Figure 18: Detail of the limiting curves B for the g-state Potts antiferromagnet on 
a square lattice 7p x oo^. (a) Region around the T point at q « 2.737 + 1.405 i 
and the gap between q « 2.6590 + 1.525 i and q « 2.7275 + 1.4173 i. (b) Detail 
of the region around the T point. At this point we have t w (0.125,0.581,0.425) 
and 9 w (0.248,1.053,0.804). (c) Region around the gap between the points q ~ 
2.7619 + 0.46936 % and q w 2.7873 + 0.47546 z. 
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Figure 19: Zeros of the partition function of the g-state Potts antiferromagnet on a 
square lattices 8p x 40p (squares), 8p x 80f (circles) and 8p x oop (solid line). The 
isolated limiting zeros are depicted by a x. The limiting curve was computed using 
the resultant method. 
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Figure 20: Detail of the limiting curves B for the g-state Potts antiferromagnet on a 
square lattice 8p x oo F . (a) Region around the gap at q ~ 1.39 + 2.59 i. (b) Region 
around the T points q « 2.783 + 1.088 i and q « 2.801 + 1.043 i. At the former T 
point we have t « (1.074, 0.0184, 1.115) and 6 w (1.642, 0.037, 1.679); at the latter we 
have t w (0.119,0.813, 1.031) and 9 w (0.237, 1.365, 1.601). (c) Detail of the T point 
at g ~ 2.783 + 1.088 i. (d) Region around the tiny gap at q ~ 2.7515. 
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Limiting curves for the square-lattice strips Lp x oof with 3 < L < 8. 
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